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ABSTRACT 


Lukens, Larry Allan, Ph.D., Purdue University, 
March 1971. An Experimental Investigation of Electric 
Field Intensity and Wall Heat Transfer for the Heating 
Region of a Constricted Arc Plasma. Major Professor: 
Miank P. Incropera. 

Electric field intensity and total wall heat 
transfer measurements have been performed for the inlet, 
asymptotic, and anode regions of a laminar, atmospheric 
cascade arc. In addition to the diagnostics required 
for these measurements, a unique radiation gage has 
been developed and used to obtain radiative wall heat 
flux data. While the majority of the measurements have 
been made for an argon arc, some data are reported for 
helium and nitrogen arcs. 

Comparisons between the experimental data and 
available equilibrium and nonequilibrium numerical 
solutions are presented in an attempt to ascertain the 
adequacy of single-fluid models to accurately describe 
Overall arc energetics. For the most part, the equili- 
brium theory predictions for the asymptotic region in an 
argon arc agree with the experimental results to within 


the combined experimental and theoretical uncertainties. 


However, there is a small, but discernable nonequilibrium 





XVL 


effect at both the lower and higher are currents in argon. 
In contrast, the nonequilibrium theory predictions are 

in excellent agreement with the experimental results over 
the entire range of asymptotic conditions investigated in 
argon. Further comparisons are presented which indicate 
that atmospheric helium arcs are characterized by a signi- 
‘ficant departure from local thermodynamic equilibrium. 

On the basis of comparisons between experiment and 
theory for the anode eiion, it is shown that anode heat 
balance models are ill-suited for providing design es- 
timates for the total heat flux to the anode of the arc 


constrictor. 











GHAETER I se INTRODUCTION 


I.l General Objectives of This Study 
Owing to interest in space flight technology, MHD 


systems, and thermonuclear energy conversion, considerable 
effort has been devoted to the study of high temperature 
gas flows. Of particular interest to this study is that 
portion of this effort which deals with a device which has 
long served as a versatile tool in generating high tempera- 
ture gases for both industrial and laboratory use - the 
plasma arc generator. 


Designed initially by Maecker 1] 


, the operation of 
the Gerdien, or constricted tube, arc generator is based 
upon the conversion of electrical to thermal energy. In 
the most common configuration (Figure I.1), this conversion 
is accomplished by passing a gas through an electric arc 
which is established by applying a d-C potential between 
the electrodes. 

The constricted tube arc produces a continuous and 
well-defined plasma flow which, in addition to its use in 
certain technological applications, is particularly useful 


for the experimental determination of high temperature 


gas transport properties and the development of high 





temperature diagnostic techniques. Irrespective of its 
use, however, the constricted tube arc is distingumenea 

by extremely high heat fluxes to the containing walls and 
these walls are often destroyed if sufficient cooling is 
not provided. Accordingly, an understanding of the elec- 
trical and thermal characteristics of arc constrictors 

is important to the designer as well as to those with a 
general interest in arc physics. The objective of this 
study is to enhance this understanding through an experi- 
Memea: investigation of the electrical and thermal effects 
in an atmospheric cascade arc and through an interpretation 
of these effects in terms of available theoretical studies. 
While the major portion of the experimental work has been 
performed for an argon plasma, some measurements have 

also been made for nitrogen and helium arcs. 

The regions of interest for the cascade arc are 
depicted schematically in Figure I.1l. The arc column, 
which is a composite of individually water-cooled and 
electrically isolated copper discs, is divided into three 
distinct regions: inlet, asymptotic, and anode. The gas 
flow in each region is subjected to an externally applied 
electric field, and the inlet and asymptotic regions are 
generally bounded by non-current-carrying surfaces. In 


contrast, the anode surface is electrically conducting. 
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Figure I.1 Designation of Arc-Heating Regions 





The inlet region is further distinguished by a significant 
axial variation of the flow variables, whereas the asymp- 
metic region is characteristic of thermally fully-devel- 
oped flow. As a matter of convention, the axial location 
at which the mean (mass averaged) enthalpy attains 90% 

of its fully-developed value is chosen to mark the be- 


ginning of the asymptotic region. 


ieee Pvevrous investigations of 


Plasma Tube Flow 
Although the literature devoted to theoretical study 

of the inlet and asymptotic regions is extensive, few 
experimental studies have been performed for purposes of 
determining the overall arc energetics. Of particular 
scarcity are experimental studies pertaining to the inlet 
region of the water-cooled duct. pieseaee performed 
electrical potential and wall heat flux measurements in 

the inlet region of a helium arc, and investigators at 


[3-5] have made similar 


the NASA-AMES Research Center 
measurements in air and nitrogen arcs. However, the range 
of experimental conditions (gas flow rate, arc current, 
and duct geometry) investigated in these studies differ 
Significantly from those of this investigation. Conse- 


quently, meaningful comparisons cannot be made between the 


data of this study and the experimental results reported 





by these investigators. Similarly, Pfender, et moe) 


report results for the local wall heat flux in the inlet 
region of an argon cascade arc; however these measurements 
were for turbulent flow and were used primarily for the 
support of anode heat transfer studies. 

A greater number of experimental investigations 
have been performed for the asymptotic region of the 
ee While these studies are primarily concerned 
with the use of an inversion method to determine plasma 
transport properties, they necessarily include measure- 
ment of the arc-column electrical characteristics and 
are therefore germane to this study*. In addition, Morris, 
et ee present asymptotic wall heat flux data and 


[15] 


Emmons [13] and Maecker report on the percentage con- 


tribution of radiation to the total energy balance. 
In view of the excellent reviews given by 


Skifstag!16! sa 


and, more recently, by Bower the present 
discussion of theoretical investigations is limited to 
those studies which are to be used for comparison with 
experiment in Chapter IV. 

The approximate analytical solution by Stine and 


[18] 


Watson was one of the first reported theoretical 


treatments of the inlet region of the cascade arc. This 





- The inversion method of property determination will be 
discussed further in Chapter IV. 





solution is based upon a reduced form of the governing 
energy equation, obtained by neglecting radiation, radial 
convection, and viscous dissipation. In addition, the 
Stine-Watson model assumes a uniform mass flux throughout 
the constrictor and uses linearized relationships between 
the enthalpy and related properties (T, k/C,s and o). The 
primary weakness of this model is its failure to properly 
account for the variation of the transport properties 
with temperature and, particularly when applied to the 
higher temperature argon and nitrogen arcs, its failure 
to account for radiative losses. 

To remove the need for the above assumptions, de- 
tailed numerical solutions for the inlet region have been 


[19] 


carried out by Watson and Pegot and by Bower and 


Incropera |29! | 


Of the two, the latter will be used for 
the comparisons of this study since it gives greater con- 
Sideration to pertinent thermal effects. Specifically, 
Bower and Incropera formulate a rigorous equilibrium model 
for the laminar, steady, and axially symmetric (without 
swirl) flow of an atmospheric argon plasma through the 
heating region of a constricted tube arc. An implicit 


[17] is then used to solve the 


finite-difference scheme 
resulting conservation equations (mass, momentum, and 


energy) consistent with the thermal and caloric equations 





of state, Ohm's law, and the temperature dependence of the 
required thermodynamic and transport properties. Implicit 
in the final form of the equations are the boundary layer 
approximations and the assumptions of negligible radial 
current density, quasi-neutrality, and local thermodynamic 
equilibrium. The results of this numerical solution are 
presented in the form of wall heat transfer and shear 
stress (as predicted from energy and momentum balances) and 
correlations between the mean Nusselt number and an Ohmic 
heating parameter. 

As a result of the evolution of the aforementioned 
numerical solutions, it would appear that the theoretical 
treatment of the inlet and asymptotic regions of arc 
constrictors has reached a practical limit. However, the 
variable tube diameter experiments performed by Emmons !133 
and, more recently, the spectroscopic measurements con- 


ducted by Bott [1! aes 


and by Incropera and Giannaris 
provide strong experimental evidence of the presence of 
nonequilibrium in confined plasmas. Accordingly, the 
validity of the local thermodynamic equilibrium assumption 


which is implicit in the single-fluid model is subject to 


question.* 





* A comprehensive discussion of local thermodynamic [22] 
equilibrium in confined plasmas is given in references 


and L23] 





In an attempt to account for the nonequipartition 
of species translational energies (thermal nonequilibrium) , 
Pytte and Winsor |?4 formulated a nonequilibrium theory 
for the asymptotic region of a helium arc. This attempt, 
however, is at best, semi-empirical, for its solution 
relies upon experimentally determined arc-column charact- 
eristics. In addition, the Pytte-Winsor model suffers 
from use of the Saha (equilibrium) equation to compute 
the plasma composition. 

While accounting for both thermal nonequipartition 


[25] has formulated a 


and chemical nonequilibrium, Clark 
rigorous multifluid (electrons, ions, and neutral atoms) 
model for the steady, laminar flow of a nonequilibrium 
argon plasma in the inlet and asymptotic regions of a 
constricted tube arc. After applying the boundary layer 
approximations, the governing set of species and global 
conservation equations, along with the thermally perfect 
gas equation of state and Ohm's law, are solved numerically 
using finite-difference techniques 2°! Of par ereular 
interest to this study are the predictions for the asymp- 
totic electric field intensity, total wall heat flux, and 
ratio of radiative to total wall heat flux which result 
from the nonequilibrium solution. 


The need for extended operating times for the most 


Severely heated component of the arc constrictor, the 





anode, has prompted considerable research on the complex 
energy exchange mechanisms occurring at this surface. 
Energy transfer models have been proposed in which the 
various contributions to the total anode heat flux are 
delineated. For the most part, these models are based on 


[Glee 2 oes 


a heat balance across the anode surface 
addition, extensive studies of the electrical characteris- 
tics of anodes have been conducted by researchers at 


iemouth ?1-33! . 


For purposes of comparison with the 
anode heat transfer measurements of this study, primary 
emphasis will be placed upon the widely used model which 


is described by Pfender, et a1. 161, 


I.3 Scope and Specific Objectives of This Study 


The primary objective of this study is the measure- 
ment of the axial distribution of the electric field 
intensity and the total wall heat flux for the inlet, 
asymptotic, and anode regions of an atmospheric argon 
cascade arc. These measurements are made for a wide 
range of laminar flow operating conditions. Additional 
Measurements are made to determine the ratio of radiative 
to total wall heat flux for the asymptotic region of an 


argon arc**., An important secondary objective includes 





* A review of the subject of anode heat transfer for 
various electrode geometries is given by Eckert and 


Pfender [39! 
** A review of the literature pertinent to this subject is 
presented in the Appendix. 





10 


comparison of the asymptotic wall heat transfer and electric 
field intensity measurements with the experimental data of 
other investigators and with the analytical and equili- 
brium and nonequilibrium numerical solutions described 
previously. With respect to the latter, particular atten- 
tion is given to ascertaining both the validity of existing 
Nusselt number correlations and the adequacy of single- 
fluid models for the description of arc energetics. Pur- 
suant to gaining additional information on anode physics 
and the suitability of anode heat transfer correlations, 
further comparisons are made with existing anode heat 
transfer models. The measurements are also of value in 
providing an improved understanding of the development 
length characteristics of the constricted tube, laminar 
arc. In order to ascertain the influence of operating 

gas, additional experimental results are presented for 
nitrogen and helium arcs. 

This study, which is the first serious attempt to 
experimentally verify existing equilibrium and nonequili- 
brium theoretical predictions of thermal parameters for 
the arc heating region, is intended to provide information 
useful to the arc designer and to contribute to an im- 


proved understanding of arc physics and flow modeling. 





Heal 


CHAPTER II ARC GENERATOR AND DIAGNOSTIC TECHNIQUES 


II.1l Arc Generator and Supporting Systems 
A well-defined, stable plasma flow is achieved 
through use of a constricted tube plasma arc generator 


[34] (Figure II.1). The 


manufactured by Creare Inc. 
constrictor is formed by stacking together a series of 
individually water-cooled copper segments (discs), each 

of which is 0.632 cm wide with a 1.0 cm central bore. Each 
segment is electrically insulated from adjoining segments 
by cool O-ring seals which further provide for alignment 
and hermetic sealing of the duct. The lengths of the 

arc heating and field-free regions may be varied according 
to the anode location and the number of segments used to 
form the stack (up to a maximum of 40). The thoriated 
tungsten cathode 1s water cooled and may be precisely 
centered with respect to the initial nozzle-shaped seg- 
ment through micrometer adjustment. The cathode-nozzle 
geometry provides for a high degree of flow symmetry which, 
along with water cooling, provides for extended cathode 
life. If the tungsten tip of the cathode is periodically 


(every 20 hours) re-shaped on a lathe such that the proper 
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cathode-nozzle geometry is maintained, a single cathode 
will provide approximately 100 hours of operation. 

The supporting systems which are necessary to sus- 
tain the arc generator are illustrated schematically in 
Figure II.2. Electrical power is supplied by a D-C 
motor-generator set capable of delivering 75 kw at 300 
amperes. Stable operating conditions are achieved through 
use of two, teewice variable ballast resistors which 
are connected in series with the power supply and arc 
generator. The electrical system will allow for the use 
of either r-f discharge or prod starting procedures. 

The main coolant system supplies cooling water to 
the plasma generator coolant manifold and to the main 
ballast resistor. The manifold in turn distributes cooling 
water to the cathode, the starting prod, those segments 
which do not serve as calorimeters, and to the additional 
downstream ducting required to cool the gas for discharge. 
A secondary coolant system is employed to reduce flow 
demands on the main coolant system and to provide a 
means of selectively monitoring and controlling the 
coolant flow rates through those duct segments for which 
wall heat transfer measurements are desired. This system 
also supplies cooling water to the auxiliary ballast 


resistor. 
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Receiving its supply from commerical gas cylinders, 
the gas system facilitates control and measurement of the 
operating gas flow rate. Stringent purity requirements 
are maintained throughout the gas system including use of 
research grade gas, prior ultrasonic cleaning of all 
components exposed to the plasma stream, and continuous 


outgassing of the system prior to operation. 


Pie) ees bewlar Operating Procedures 


and Limitations 

Despite the aforementioned purity requirements, the 
level of duct contamination resulting from arc operation 
has steadily increased over the past two year period. 
During the earlier stages of this research, when the major 
portims of the electric field intensity and total wall 
heat flux data were acquired, there were no visible indi- 
cations of surface contamination, even after extended per- 
i1ods of operation (up to 40 hours). However, during the 
latter phases of the research, when the radiative wall 
heat flux data were acquired, the level of contamination 
continually increased in spite of special precautions that 
were taken in addition to those listed above. 

The contamination is evidenced by a thin deposit 
which coats the duct walls in a distinct and repetitive 


pattern. The deposit, which appears to be a dark oxide 
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layer, forms on the leading (upstream) edge of each of the 
segments in the arc heating region and diminishes in the 
axial direction, leaving the downstream half of the seg- 
ment surfaces free of visible contamination. The deposit, 
which at times has a bluish cast, has not been observed in 
the field-free region, not even on the segment immediately 
downstream of the anode. 

Similar deposits have been observed by other in- 


i398 74.33, 


vestigators however opinions differ concerning 


the effect of surface contamination on experimental results. 


el 3 


Emmons claims that the overall arc characteristics are 


virtually unaffected by small amounts of contamination. 


On the other hand, meee 32! 


reports a reduction in 
the electric field intensity for helium of some 10% when 
data are acquired under contamination-free conditions. 
Throughout the course of this study, both electric field 
intensity and total wall heat flux measurements have been 
consistently reproducible to within their respective 
experimental uncertainties, irrespective of the level of 
contamination. None the less, it was assumed initially 
that contamination would render the proposed radiation 
gage completely useless, since the dark deposits would 


Significantly alter the optical properties of the gage 


surface. While useful data were acquired from the 
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radiation gage, the contamination level within the duct 
still prevented the gage from being used as it was 
originally intended (see the Appendix). 

As reported by McKee, et a1, 132), the difficulties 
associated with obtaining reproducible measurements of 
anode characteristics are most likely due to the pro- 
nounced affect that surface contamination has on electrode 
performance. In connection with this problem, a number 


beiieea oN have devoted considerable 


of investigators 
effort to the elimination and investigation of surface 
contamination within constricted tube arcs. In summarizing 
work performed in this area at Dartmouth, pidecagwes 
Suggests a number of procedures which may be used to re- 
duce the level of contamination within the duct. These 
recommended procedures include use of ultra-pure grade 
gas, elimination of possible cathode contaminants through 
use of a sonic orifice, outgassing of the duct surfaces, 
and nickel plating of all copper surfaces exposed to the 
plasma stream. The plating procedure was found necessary 
in order to eliminate outgassing from the copper surfaces 
during arc operation. 

Similar, yet considerably less involved, measures 
were taken in this study in an attempt to eliminate the 


contamination of the duct surfaces. Ultra-pure grade 


(as opposed to commercial grade) gas was used; however 
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meense did not result in a significant reduction in 
surface contamination. This is not to say, however, that 
use of commercial grade gas will not result in surface 
contamination. After a thorough cleaning of the gas 
system failed to effect any change in the degree of sur- 
face contamination, the inlet portion of the gas system 
was virtually eliminated by connecting the gas cylinder 
directly to the plasma generator (a short piece of plastic 
tubing was connected between the gas inlet to the plasma 
generator and the gas regulator and the regulator was 
connected directly to the gas cylinder). In addition, 
the gas regulator (a two-stage, metal diaphram type) was 
completely disassembled, cleaned ultra-sonically with 
acetone and isopropyl alcohol, and allowed to degas at 


ion ° 


mm Hg for a period of 24 hrs. Several hours prior 
to starting the plasma generator, the combined generator 
and gas system was evacuated to approximately 40 microns 
and the metal components, particularly the duct segments, 
were heated to accelerate degassing. Unfortunately, none 
of, or any combination of, these procedures led to any 
Significant reduction in the level of surface contamina- 
tion. In addition, the possible influence of prod and 
r-f discharge starting techniques was investigated, as 


were the procedures used to mechanically and chemically 


clean the duct segments. Special care was taken in 
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positioning the cathode within the initial nozzle segment 
so as to minimize possible contamination from erosion of 
the tungsten cathode tip. Again, there were no per- 
ceptable changes in the level of surface contamination with 
any of these additional measures. 

On the basis of these results, it is tentatively 
concluded that the surface contamination is an oxide 
layer which is formed by a chemical reaction between the 
copper duct walls and the gases which escape from these 
walls during operation of the arc. This conclusion is 


2 concerning 


supported by the observations of Runstadler 
the nickel plating of the copper surfaces, and also by 
the experiments performed in this study with the thin 
film aluminum surface on the radiation gage (Appendix). 
Further support of this conclusion is provided by the 
x-ray diffraction analysis conducted by McKee, et ai, '32! 
which identified the surface contaminants as copper-oxide 
and carbon. While not fully understood, the form of the 
axial distribution of the surface contamination on seg- 
ments in the arc heating region and the absence of surface 
contamination in the field-free region indicate that 


certain electrical effects must also influence the 


contamination. 
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In addition to experiencing contamination problems, 
some difficulties arose in connection with generating He 
and N, arc plasmas. The maximum current-carrying capa- 
bilities of the starting prod (75 amp) and r-f discharge 
(50 amp) circuits, while adequate for initiating argon 
arcs, are insufficient to permit starting in either 
nitrogen or helium. Accordingly, a special pre-start 
procedure was devised. This procedure involves use of a 
dual gas supply manifold whereby both Ar and either No or 
He can be supplied to the arc generator simultaneously. 
Following initiation of a stable, 50 amp Ar arc, the 


proportion of N. or He to Ar is slowly increased from 0 to 


Z 
100% while the open-circuit voltage is steadily increased 
to Maintain a stable operating condition. While this 
procedure is adequate for initiating No and He arcs, the 
maximum obtainable open-circuit voltage of the power 
supply is nonetheless insufficient to sustain stable 
Operation of these gases over as wide a range as is 
possible with Ar. 

The restrictions which are placed upon both high 
and low current operation of the arc generator are a re- 
sult of certain physical limitations of the power supply 
and external ballast resistor. Specifically, high current 
Operation is limited by the maximum open-circuit voltage 


Capability of the power supply. In contrast, low current 


Operation is restricted by the lower limit of electrical 
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@eability which, in turn, is a function of the available 
ballast resistance as well as the maximum obtainable open- 
circuit voltage of the power supply. 

These limits of operation can be explained by con- 
sidering a set of hypothetical (yet typical) load line 
characteristics for an arc generator and its associated 
D-C power supply (Figure II.3). As indicated by the 
electrical circuit which accompanies Figure II.3, the 
power supply load line characteristic is a plot of the 
voltage available to the arc generator versus the available 
arc current. Generally, the power supply characteristic is 
represented by a straight line which intersects the vol- 
tage axis at the value of the open-circuit voltage and 
which has a slope equal to the negative of the sum of 
the constant internal resistance of the power supply (Ro) 
and the variable external ballast resistance (R)) - The 
arc generator load line, on the other hand, represents a 
plot of the current conducted by the arc versus the 
voltage required to sustain that current. On the whole, 
the arc generator is characterized by a wide spectrum of 
load lines which result from various combinations of the 
Operating gas, the operating gas flow rate (m), the 
Operating pressure (p), the overall (cathode to anode) duct 


length (L), and the duct radius (R). 
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The actual operating state (for a specified gas, 
hm, p, L, and I) of a constant diameter arc generator is 
uniquely determined by the point of intersection of the 
load (arc generator) and supply characteristics. The 
degree of electrical stability which exists for a given 
@eeracing state is a function of the relative slopes of 
the two load lines at their point of intersection. The 
criterion that must be satisfied in order to achieve 


stable operation is given by 


IV 


0 (Tiga) 
LOAD SUPPLY 


which, in the case of a linear supply characteristic, 


reduces to 


(==) eee 0 (Eie 2) 


It is noted at this point that, while Equation II1.2 
establishes, in principle, the criterion for stable opera- 
tion, in practice the slope of the load characteristic must 
always be greater than the slope of the supply character- 


istic if stable operation is to be realized. That is, 
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(=) aR + Rage > 0: (ieee: 3) 
~ LOAD 
' : aE ; 
Ernce Re 1S constant and since (a> 1s dictated by the 


LOAD 
desired arc operating condition, improved electrical sta- 


bility at each condition can be effected only by increasing 
R,- This, in turn, necessitates an increase in the open- 
circuit voltage. Accordingly, for stable arc operation, 

it is generally advisable to use the highest possible 
settings of both the ballast resistor and the open-circuit 
voltage of the power supply. 

Irrespective of the value of Rye Beat ioneelt . 3 cwaull 
be satisfied as long as operation is confined to the 
positive slope portions of the arc generator load iemes 
Nevertheless, the maximum open-circuit voltage capability 
of the power supply can place a restriction on high cur- 
rent Operations. For example, if all of the ballast 
resistance has been removed in order to operate in condi- 
tion-a (Figure II.3), further operation at higher arc 
currents (position-a') or operation at greater values of 


m and/or L (position-a") can be achieved only by increasing 





* Note that operation at the initial point of intersection 
(e.g., points x and y in Figure II.3) is not likely since, 
according to Equation II.2, the arc would be inherently 
unstable at this condition. 
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the open-circuit voltage. Therefore, aside from possible 
over heating of the electrodes, high current operation 

is limited by the maximum obtainable open-circuit voltage 
of the power supply. 

According to the stability requirement expressed by. 
Equation I1I.3, operation on the negative slope portions of 
the arc generator load lines (e.g., position-b) is in- 
herently less stable than operation in the positive slope 
regions. In addition, a lower limit for stable operation 
exists at the current for which the slopes of the load and 
supply characteristics are equal (e.g., position-c). In 
practice, the arc will extinguish prior to this lower 
limit; however, lower current operation (position-c') or 
Operation at higher valves of th and/or L (position-c") can 
be achieved by increasing the open-circuit voltage and 


increasing R This procedure can be continued until 


b* 
the maximum open-circuit voltage of the power supply (or 
the maximum value of Ry) is reached. 

For the duct lengths and operating gas flow rates 
employed in this study, the point at which the slope of 


the load characteristic is equal to zero cee = 0) 
di 
LOAD 
occurs at an arc current of approximately 25 amp for Ar, 


150 amp for N and at a current in excess of 400 amp in 


oY 


He. Therefore, according to the previous arguments, low 


current operations in all three gases and high current 
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operations in N, and He are limited by the power supply. 
High current operation in Ar is limited by the maximum 
power that can be dissipated by the main ballast resistor. 
It is further noted that over heating of the anode (as 
evidenced by ablation of the anode surface) occurred at 


the higher arc currents in’N, (I > 200 amp) and in He 


Z 
=> 100 amp). 


in. 3 Diagnostic Methods and Related Instrumentation 


The overall and local arc-column characteristics 
are determined by measuring the potential difference 
between the "floating" duct segments and the cathode. The 
arc voltage measurements are used both to determine the 
axial distribution of the electric field intensity and, 
along with the arc current and calorimetric measurements, 
to determine the mixed-mean enthalpy of the gas at each 
segment location. The arc voltages are measured with a 
high input impedance (4 megaohm) strip-chart recorder 
(Electronik-19 by Honeywell). The arc current is deter- 
mined by monitoring the voltage drop across a calibrated 
shunt resistor with a digital voltmeter (Dana, model 
# 54030). As shown schematically in Figure II.4, further 
circuitry is available to permit up to 400 ma of current 
to be drawn from a duct segment while its potential (with 


respect to the cathode) is monitored. 




















Figure II.4 
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The energy loss to the duct segments, to the 
cathode, and to the remaining arc generator components is 
determined calorimetrically. Accordingly, each of these 
components is equipped with a specially designed tempera- 
ture difference (AT). transducer. The AT transducers 
(Figure II.5) employ a single pair of copper-constantan 
thermocouples, mounted in plexiglass inlet and outlet 
plenum chambers, and wired in series to provide a voltage 
Signad which is directly proportional to the coolant 
temperature rise. Considerable care was taken in the 
Shielding and calibration (against mercury-in-glass 
thermometers) of these transducers and they are accurate 
to within 5% of the indicated reading. The transducer 
output voltages are amplified (Dana model #3850-V2) prior 
to being displayed on the digital voltmeter. The abso- 
lute temperature of the operating gas at both inlet and 
discharge locations is measured with a copper-constantan 
thermocouple (as referenced to a copper-constantan ice- 
junction). The gas temperatures are monitored by the 
digital voltmeter. 

A unique radiative heat transfer gage is used to 
delineate the radiative contribution to the total wall 
heat flux. As shown in Figure II.6, the radiation gage 
consists of two semicylindrical copper segments which 


are thermally isolated from one another, individually 
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water cooled, and supplied with individual AT-transducers. 
The inside wall of each half segment is coated with either 
a highly reflecting (aluminum) or highly absorbing (black) 
surface. The overall dimensions of the assembled gage 

are precisely the same as those of a standard duct seg- 
ment. A thorough discussion of the principle of opera- 
tion, the fabrication, and the associated experimental 
problems involved in the use of the radiation gage are 
left to the Appendix. 

The operating gas flow rate and the coolant flow- 
rates through the cathode and those duct segments serving 
as calorimeters are measured with precision accuracy 
(+ 1%) flowmeters (Fischer & Porter model no. 10A0735). 
The two gas flowmeters have a combined range of from 0.02 
to 1.50 1lb/min in argon. At the time of each flowmeter 
reading, the inlet gas temperature and pressure are 
measured and used to correct the indicated flow rate 
should they differ from the temperature and pressure of 
calibration. An additional precision accuracy flowmeter 
is used to measure the total coolant flow rate to the 
plasma generator. 

With measurement of the arc current (I) and axial 


voltage distribution (E), the gas inlet temperature (T ) 


INLET 


and flow rate (m), and the power losses to the duct 


segments (PLS) and to the cathode (PLC), the axial 
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distribution of the gas enthalpy can be determined by 


performing energy balances at successive axial locations. 


Z(PLS) +£ 
ry - 
<2 


Mm, Titer 


PLC +4: 





Figure II.7 Energy Balance Diagram 


With reference to Figure II.7, the enthalpy of the gas 
at the exit plane of a segment in the arc heating region 


1s given by 
: 
i aoe (PLS) . - Ee 
J. i=l 


h. ms eS Tinrer * 2 : (Eig) 
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This value is then algebraically averaged with the enthalpy 
determined by the previous (j-l) energy balance, and 

the resulting "mixed-mean" enthalpy is assigned to an 

axial location corresponding to the center of the 4th 
segment. 

To experimentally verify the accuracy of the diag- 
nostic procedures employed in this study, an attempt was 
made to "close off" the energy balance for the entire 
plasma generator. While the precision of this energy 
balance is limited by unavoidable ambient losses, it was 
still possible to obtain closure to within the cumulative 
experimental uncertainty of the component measurements 
(t 20%). A more precise check is one based on a local 
energy balance performed on a segment located in the 
fully-developed portion of the arc heating rection: « In 
every case, closure of this local energy balance was 
obtained to within 5% - a value comparable to the total 
uncertainty of these measurements (+5-10%). The consistent 
reproducibility of data taken over a two year period pro- 
vides further confidence in the accuracy of the measure- 
ments as well as the consistency of the overall system. 


A detailed uncertainty analysis is used for estima- 


ting the uncertainties associated with each measured 





* The nature of this energy balance will be discussed 
in detail in the following chapter (Section III.3). 
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variable and for calculating the propagation of these 
uncertainties into each reported result. The method used 
in this analysis is the standard approach for single-sample 


ee The uncertainty calculations, as well 


experiments 
as the conversion and reduction of raw data, are per- 
formed on the computer. The thermo-physical properties 


required for determining the mean temperature, Tan mm) 


and the Nusselt number-Ohmic heating parameter correlations 


_ and a) are obtained from seueto!>”! and Drellishak, 
Knopp, and eenpoe in the case of Ar; Lick and pumens >?! 
and DeVoto and pi 40) in the case of He; and Drellishak, 
Aeschliman, and Camber /41 and Aveo! for N.. The re- 


2 


ported uncertainties of these properties are also consi- 
dered in the aforementioned uncertainty analysis. A more 
detailed description of the arc facility, diagnostic tech- 


Nigues, and uncertainty analysis is given by aense >? . 
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CHAPTER III EXPERIMENTAL RESULTS 


Telia): Range of Measurements 


Experimental data were acquired for Ar, No, and 


He arcs over the following ranges of laminar flow operating 


conditions: 
Argon; I = 35 - 240 amp, m = 0.03 - 0.11 1b/min, 
constrictor length = 5 in (20 segments). 
Nitrogen; I = 100 - 200 amp, mh = 0.03 lb/min, 
constrictor length = 5 in (20 segments). 
Helium; I = 50 amp, m = 0.07 lb/min, 


constrictor length = 2.5 in (10 segments). 


All measurements were made with a 1.0 diameter copper 
duct operating at near atmospheric pressure. The assertion 
of laminar flow is based upon the flow transition studies 


oe and by Runstadler '~"! in 


by Incropera and Leppert 
which the onset of turbulence is reported to occur ata 
mean Reynolds number (based on properties evaluated at T 
of 700 to 900. The largest value of the mean Reynolds 


number encountered in this study is approximately 600 


(Argon, I = 35 amp, m = 0.11 lb/min). Hence, the above 
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stated ranges of operating conditions are assumed to be 
representative of laminar flow. 

Rather than presenting a complete tabulation of the 
data acquired for the inlet region, certain data were 
selected to demonstrate the effect of varying gas flow 
rate and arc current on the development characteristics 
of the arc. All of the acquired data pertaining to the 
asymptotic region are presented along with related data 
from other investigators. The interpretation of the 3 
experimental results as well as comparisons with related 


theories will be the subject of Chapter IV. 


III.2 Electrical Characteristics 
Prior to presenting the experimental data for the 
arc-column (electrical) characteristics, a discussion is 
in order concerning the method employed in the acquisition 
of this data. 
A method frequently used in the measurement of 
arc-column characteristics involves use of the duct seg- 


[10,12,15,32] rt 


ments as "floating" potential probes 
1s commonly recognized, however, that a floating duct wall 
will assume a potential which is negative with respect 


to the plasma core by an amount equivalent to the plasma 


* 
sheath potential . Accordingly, on the assumption that 





* Further treatment on the formulation and characteristics 
Of plasma sheaths is given in references [30] and [45]. 
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the sheath potential is appreciable, some investigators 
have questioned the suitability of using a floating duct 
wall to measure plasma potentials. 

In reference to their cascade arc potential measure- 
ments, this question has been raised most strongly by 


jee? In performing their measurements, 


Pfender, eta 
they allowed a current (I) to be drawn from the duct 

wall while they simultaneously monitored the potential of 
this wall with respect to the cathode. These measure- 
ments indicated that the recorded wall potential was 
strongly influenced by the magnitude of the current drawn. 
In particular, they observed that the voltage increases 

by 40 to 100% over the floating wall (I. = 0) potential 
when a current of up to 50 ma is drawn. Above 50 ma, the 
segment potentials were observed to be independent of the 
Current drawn. With regard to this observation, they 
argued that, for Ir. > 50 ma, the charge separation which 
otherwise would exists in the plasma sheath region is 
exactly cancelled and, consequently, the duct wall is at 
the potential of the plasma column. They therefore con- 
Cluded that floating wall potential measurements do not 
provide a suitable indication of the actual plasma 
potential, but enee accurate measurements of this para- 
meter can be made by recording the wall potential while 


a current in excess of 50 ma is drawn from the wall. 
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[12,13] indicates that the sheath 


In contrast, Emmons 
potentials for an argon-copper duct combination are 
approximately 1 to 2 volts, for which case the floating 
wall potentials are suitably close to those in the arc 
column. In support of this contention is the experimental 


ed 


technique used by Maecker wherein every attempt was made 


to eliminate impressed wall currents a) while performing 


146] on the 


plasma potential measurements. Schreiber 
other hand, has measured significant voltage increases 
With current drawn, but he attributes these increases to 
voltage drops across the relatively cool gas boundary lay- 
ers. Additional support for Schreiber's deduction can 

be inferred from the measurements made by researchers at 


[32733] While the Dartmouth data were used 


Dartmouth 
primarily in support of anode studies, the arc configura- 
tion and external circuitry used for these measurements 
are similar to those employed in this study (Figure II.4) 
and to those used by Pfender, et aie 

In an attempt to resolve these inconsistencies, the 
Variation in the duct wall potential with current drawn 


was measured at various axial locations and arc currents. 


As shown in Figure III.1, these measurements indicate that 





* Maecker contends that the measured wall potential 
difference between two segments will be in error by an 


amount proportional to the impressed current, I: 
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the percentage increase in wall potential with increasing 
I diminishes both with increasing axial distance and 
with increasing arc current. They also indicate that 

the wall potential increases continuously with increasing 
I, over the entire range of T values investigated. In 
this respect, the distributions of Figure III.1 are 
qualitatively similar to the potential distributions 


Lt] and to the lower current ranges 


[33] 


recorded by Maecker 
of anode conduction reported by Runstadler 
Simultaneous calorimetric and wall potential mea- 
Surements indicate that the total wall heat flux to a 
duct segment does not vary appreciably (< 3%) over the 
entire range of I.: In addition, the value of the 
asymptotic electric field intensity (Ba! was observed 


to be independent of the current drawn. This fact 


joi) 1, [32] 


is discussed by Emmons and Land and by McKee, et a 
and is strongly supported by the experimental measure- 
ments of Morris, et ai, 241, wherein the voltage gradient 
was measured using three independent methods. With the 


first method, the value of EY was computed from the 


a 
/ 
asymptotic energy balance equation (this equation is dis- 
cussed later in this chapter) along with calorimetric 
Measurements of the power losses to the duct segments. 


The second method, which is similar to that employed by 


Maecker, utilizes two of the duct segments as potential 
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probes. Unlike Maecker's method, however, a small current 

was supplied to the arc through these probes (in an attempt 
to cancel the voltage drop across the "contact resistance" 

between the probes and the arc column) while the potential 

difference between the two duct segments was measured. 

This potential drop was then divided by the distance be- 


tween the two probe segments to give the value of E The 


xja- 
third method used by Morris, et al. involved a cathode 
withdrawal experiment wherein the total arc current was 
held constant while the total (cathode to anode) voltage 
drop was recorded as a function of the total (cathode 

to anode) duct length. These three methods, when applied 
to the central region of a hydrogen arc, resulted in values 
of the electric field intensity which agreed to within 
approximately 10% - a value comparable to the combined 
uncertainties of these measurements. 

The agreement demonstrated by these experiments 1s 
presumably a result of the conditions which characterize 
the fully-developed region. That is, regardless of the 
magnitude of possible radial electric field components, 
the voltage difference between any two axial points in 
the fully-developed region is invariant as long as the two 
points are at the same radial location. If this were not 
the case, the axial component of the electric field 


intensity would vary with axial location and the flow 
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would no longer be asymptotic. The independence of . 
with I, was tested under the most severe condition 
(I = 50 amps, I, = 400 ma) of this investigation and, as 
shown in Figure III.2, the slopes of the linear portions 
of the wall potential distributions Cea. are the same 
fOr = = QO and I= 400 ma. It is further noted that the 
electric field intensity distribution on either side of a 
segment from which current is being drawn is unaffected 
by the wall potential variations of that segment. 

Although it is concluded that the current drawn 
from a duct segment (up to IT. = 400 ma) does not apprecia- 


bly effect the total wall heat flux or EY the increase 


a 
in wall potential with increasing T, will nonetheless have 
a direct effect on the energy balance measurements used 

to determine the mean enthalpy of the gas, particularly 

at the lower arc currents. As indicated in Figure III.3, 
the mean enthalpy of a 50 amp, argon plasma is considerably 
greater when computed using wall potentials measured with 
I. = 400 ma than it is when computed with floating wall 
potentials. While the relative accuracy of these two 
distributions could best be determined by comparison with 
experimental enthalpy profile data, such data are not 
available. However, a comparison is made with the non- 


[25] 


equilibrium numerical solution by Clark Pies 


solution, as applied to the asymptotic region of an 
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atmospheric argon plasma in a 1.0 cm diameter duct, pre- 
dicts electric field intensity and total wall heat flux 
values which agree with measured values to within the 
experimental uncertainties (< Saree As indicated in 
Figure III.3, this solution also predicts an asymptotic 
mean enthalpy which agrees well with the experimental 
value determined from the use of floating wall potentials. 
While the conclusions that can be drawn from this compari- 
son are, at best, tenuous, the agreement between theory 
and experiment shown in Figure III.3 indicates that the 
floating duct wall serves as the more suitable potential 
probe. 

After modifying the wall boundary conditions to 
include current drawn from the duct wall, the nonequili- 


aes (which includes a plasma 


brium solution by Clark 
Sheath model) was used to investigate the effect of rT. 
on the plasma sheath potential. At I = 50 amp and 

I. = 400 ma, the nonequilibrium solution predicts a 
decrease (from its floating wall value) in the plasma 
sheath potential from 3.60 volts to 1.75 volts. This 


predicted reduction in sheath potential is considerably 


less than the increases in wall potential measured under 





* Comparisons between experiment and the nonequilibrium 
model are presented in Chapter IV. 
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the same conditions. 

On the basis of the measured wall potential dis- 
tributions of Figure III.1, and on the basis of the 
Comparisons with the nonequilibrium solution, it is con- 
cluded that the observed increases in wall potential are 
primarily a result of the wall current (I) being con- 
ducted across the thin layers of relatively high electri- 
cal resistance gas adjacent to the duct ae. Accord- 
ingly, it is felt that the "floating" wall potential 
provides a more accurate measure of the actual plasma 
column potential, at least for those conditions where the 
plasma sheath potential is small in comparison to the 
wall potential (on the basis of the nonequilibrium solu- 
tion, the sheath potential, which increases with in- 
creasing arc current, attains a maximum value of approx- 
imately 5.0 volts at I = 240 amp). Therefore, all of 
the following data pertaining to the electrical character- 
istics of the arc-column are based upon floating wall 
potential measurements. 

The axial distribution of the wall potential for 


an argon arc and its variation with arc current and with 


* A decrease in the plasma sheath potential results in an 
increase in the duct wall potential measured with 
respect to the cathode. 


**This process is treated in greater detail by 


Runstadier !°7! , 
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operating gas flow rate are shown in Figures III.4 and 
TII.5, respectively. Aside from the anode fall region, the 
wall potential varies linearly with axial distance for 
entry lengths greater than 2.5 inches. The 0.25 in 

(1 segment width) location of the zero potential intercept 
is due to the positioning of the cathode within the initial 
nozzle-shaped segment. 

As indicated in Figure III.5, the increase in 
potential with increasing flow rate is to be expected 
Since, with an increase in the thermal development length 
(and thus an increase in the overall arc resistance), a 
larger potential must be applied to maintain a constant 
arc current. In contrast, if the fully-developed 
temperature profiles are, in fact, independent of gas flow 
rate (as they must be in the presence of negligible viscous 
dissipation), then the asymptotic value of the electric 
field intensity (ES a? should not vary with gas flow 
rate. Such is the case for the lower flow rates shown 
in Figure III.5. The slight variation in es with higher 
flow rates is attributed to possible departures from 
axisymmetric flow and nonuniform arc attachment at the 
anode surface. 

The asymptotic values of the electric field inten- 
mercy for Ar, No, and He and their dependence on arc 
current are plotted in Figure III.6. Included for compari- 


sOn are the experimental distributions of Runeesdier | 
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[11] Sa 


and Emmons and Land in the case of Ar, Emmons 


1, [32] [13] 


in the case of He, and Emmons 
ErS3 


and McKee, eta 


for N.. While Emmons contends that there are no 


2 
regions (in a tube with L/D = 14) for which his N. data 


2 

are independent of laminar flow rate, asymptotic condi- 
tions (as evidenced by constant values of EY and g in the 
axial direction) were observed in this study for a 
L/D z 12.6. Consequently, Emmons' lowest flow rate 
data were chosen as being the most representative of 
fully-developed flow. 

Arc-column electrical characteristics are commonly 
presented in a reduced form wherein the dependence on duct 
radius (R) is removed by plotting I/R or EI as a function 


et This reduced form is only applicable, 


or & R 
x 
however, under those conditions for which both non- 
equilibrium effects and radiation are negligible. Accord- 
ing to the variable tube radius experiments reported by 


a io and according to com- 


Emmons and Emmons and Land 
parisons made in this study between experiment and 
equilibrium and nonequilibrium numerical solutions 
(Chapter IV), conditions exist for which radiation and 
departures from equilibrium have a significant effect on 
Overall arc characteristics. For this reason, the experi- 


mental results which were selected for comparison in 


Figure III.6 pertain only to 1.0 cm diameter ducts. 
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iii swat) Heat Flux 

Typical axial distributions of the total wall heat 
flux (q) for several arc currents and laminar gas flow 
rates in argon are shown in Figures III.7 and III.8, 
respectively. As indicated in Figure III.8, gq approaches 
a constant value which is independent of gas flow rate. 
This behavior provides additional experimental verifica- 
tion in support of the assumption of negligible viscous 
dissipation. 

The asymptotic values of the total wall heat flux 
(q.) mor AY, No and He and their dependence on arc 
current are plotted in Figure III.9. In comparison with 
the values for argon, the larger values of qa (and 
a) measured is nitrogen and helium reflect the dis- 
sociation (N.) and relatively high ionization (He) 
energies which characterize these two gases. Further 
comparisons between these three gases are deferred to 
Chapter IV. 

Because of the absence of data in the literature, 
comparison of the total wall heat flux data reported 
in this study with that of other investigators is not 
possible; however, confidence in the validity of the qd. 
data in Figure III.9 may be established from the favorable 


comparisons of Figure III.6 (E, , versus I). An energy 
f 
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balance performed on a segment in the asymptotic region 
dictates that the rate of energy removal through the 
tube wall must be balanced by the local input rate of 


electrical energy. That is, 


gq. ° 27TR AX = E eee eT. (i ieee) 
a a 


According to this expression the asymptotic value of the 


total wall heat flux is expressed by, 


Ciein2 ) 





Therefore, confirmation of the qd. data shown in Figure 
III.9 is implicit with the favorable comparisons between 
the asymptotic electric field intensity data for argon 
shown in Figure III.6. Moreover, the local energy balance 
performed on a segment in the fully-developed region pro- 
vides a convenient check on the accuracy of the calormetric 
and wall potential measurements. As is evident from the 


gq. and Ey qi/2mR distributions for argon shown in Figure 
} : 


a 
III.9, closure of this local energy balance is obtained to 
within the combined experimental uncertainties of q, (+ 6%), 


E (2) 3%) , vane rte) . 
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The percentage contribution to the total wall heat 
flux due to radiation from a fully-developed argon plasma 
and its variation with arc current is shown in Figure 


IiII.10. Also shown are the experimental data of Barzelay!*?! 


ee | 


and Emmons The Harvard data are actually somewhat 


redundant in that the distribution reported by Emmons '+3! 
was obtained by applying a correction factor of 1.5 to the 
data reported by Barzelay (the total radiation meter used 
at Harvard is discussed in the Appendix). In computing 
the electrical conductivity and radiation source strength 
properties for argon, Emmons determined that the correction 
factor was necessary in order to obtain the best fit to 
overall arc data. According to Figure III.10, however, 
better agreement exists between the radiation data of this 
study and the original Harvard data reported by Barzelay. 
As discussed in Chapter II, the axial distribution 
of the mixed-mean enthalpy (h,) 1s obtained by closing off 
the energy balance at successive axial locations. The 
corresponding mean temperatures Ce are then obtained 
from the caloric equation of state [38,39,41]. Figures 
III.11 and III.12 depict the respective arc current and 
gas flow rate dependency of the axial distribution of the 
mean temperature in argon. The convergence of the data 
in Figure III.12 provides further evidence of the attain- 


ment of a fully-developed condition which is independent 


Of operating gas flow rate. 
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The variations with arc current of the asymptotic 
mean temperatures in Ar, No, and He are shown in Figure 
III.13. The extraordinarily low value of the mean tempera- 
ture in helium (nearly an order of magnitude lower than 
the corresponding value for argon) will be discussed in 
Chapter IV. 

The Nusselt number (Nu = hD/k) is introduced in an 
attempt to correlate, in terms of appropriate dimension- 
less parameters, the wall heat transfer in the asymptotic 
region of the laminar, constricted tube arc. By defining 
the heat transfer conductance (h) in terms of the total 
wall heat flux and a temperature potential, the Nusselt 


number takes the form 


2R qa 


Nu, = OT = (TT ESS 
m’m,a “w 


where KA is the thermal conductivity evaluated at Tee 
a 


and Es is the duct wall temperature (assumed constant at 
500°K)* Introduction of the local energy balance in the 
aeymptotic region (Equation III.2) into Equation III.3 


yields 


* 
While this value of T. was chosen so as to be consistent 


with the assumed value used in the equilibrium solution 
by Bower and Incropera [20], it is noted that a consider- 
able variation (say, +50%) in T. will have little effect 


on the difference (T —-T ) since T 1S Of the Ordes of 
10°*°xK. m,a WwW m,a 
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Ey a 


"aT (ih tae) 


ee 
m°om,a 


The electric field intensity is in turn related to the arc 


current by an expression of the form 


omar. 


te 
i 
NO 
3 
e9| 
O*— 


If the electrical conductivity (0) is assumed to be con- 


stant and is evaluated at soy ae then 


fEe=Si1k EF Car. Ciies 5) 
Xy 


Substituting this relationship into Equation III.4 gives 


an expression of the form 


1° 
Nu_ = ao a. (LEG) 


(TR) lade. _ mal os) 


It is now convenient to introduce an Ohmic heating param- 


eter defined as 


R CK, (5 aoed) 
a Mm m,a Ww (Teaes.7 ) 
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which may then be used to correlate the Nusselt number 
through the relationship 


me -1.0 
Nu, = 2 (Oh_) e ees) 


This correlation is plotted in Figure III.14 for argon 

along wth the correlation obtained by Bower and Incropera | *9! 
from numerical considerations. While the numerically 
predicted correlation compares reasonably well with the 
experimental data, it consistently overpredicts Nu, 
Better agreement with experiment is gained by multiplying 
Equation III.8 by a constant. In so doing the experimental 
data (for argon with I 2 100 amp) is correlated to within 
4% of the relationship 


_ 1.145 eo 
Nu, = an (Oh, ) e (TI LO) 


T 
Since Equation III.8 is independent of the type of gas, 
the above correlation should apply equally well for nitro- 
gen and helium. While the experimental values of the 
Nusselt number for these two gases (as well as for the 
lower current values in argon) fall within the range of 
values included in Figure III.14, the corresponding values 


3 > 


of Oh, fall well below (~10 ~ to 10 ~) those appearing in 


the figure. The cause of these rather marked deviations 
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from Equation III.8 is most likely due to the inaccuracies 
which result from evaluating ce and Or at the mean temper- 
ature and/or the inaccuracies which result from use of 
these equilibrium properties under conditions which are 
most likely characterized by significant departures from 
local thermodynamic equilibrium. The latter point will be 


the subject of further discussion in the following chapter. 


III.4 Anode Heat Flux 
Figure III.15 illustrates the dependency of the 
total anode heat flux on arc current for Ar, Nos and He. 
At gas flow rate and are current combinations for which 
the anode is preceded by a fully-developed region, the 
anode heat transfer was observed to vary only slightly 
e 10%) with gas flow rate and to vary linearly with are 


current (for I > 100 amp) according to the expression 


Sees (Lit «L0) 


tO 
ll 


for argon, and 


- 4300 + 367 I (Teletae | a) 


ZOr nitrogen. 





(amp) 


Arc Gurrent 
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Figure III.15 Variation of Anode Heat Transfer with 
Arc Current 
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CHAPTER IV INTERPRETATION OF EXPERIMENTAL RESULTS 


AND COMPARISON WITH THEORETICAL STUDIES 


IV.1 Discussion 


The performance of any comparison between experi- 
mental and theoretical results, of course, necessitates 
compatibility between the experimental configuration and 
the theoretical model. In this study, this compatibility 
1s achieved only for the asymptotic region of the arc. 
Although considerable data have been acquired for the in- 
let region of the constrictor, it was determined that mea- 
Surements made in this region were extremely sensitive to 
the entrance configuration. In particular, variation of 
the duct inlet geometry (effected through adjustment of 
the cathode-nozzle gap) significantly influenced the mea- 
sured axial distributions of the electric field intensity 
and total wall heat flux in the inlet region. Similarly, 
all theoretical considerations for the arc-heating region 
require the assumption of entrance flow profiles and, as 


Shown by Bower and Incropera!?9! 


, predictions for the in- 
let region depend upon the assumed form of the initial 
velocity and enthalpy profiles. Accordingly, this sensi- 


tivity of experimental and theoretical results to entrance 





ae 


conditions precludes any meaningful comparison for the in- 
let region. In contrast, the experiments of this study 


[20] indicate that entrance ef- 


and the numerical solution 
fects have no bearing on either measured or predicted quan- 
tities for the asymptotic region and, therefore, a meaning- 
ful comparison of results applicable to this region is 
possible. 

This chapter will deal only with comparisons between 
experiment and theory which apply to the asymptotic (and 
anode) regions of the constricted tube arc. In particular, 
comparisons for argon will be made between the experimental 
results of this study and the analytical solution by Stine 


[18] 


and Watson , the numerical eguilibrium solution by 


Bower and Incropera |?! 
[25] 


, and the numerical nonequilibrium 


solution by Clark Owing to the scarcity of experi- 


mental data acquired in this study for helium, the experi- 


sh will be used for comparison with 


[56] 


mental data of Emmons 


the theoretical results obtained by Bower , Pytte and 


Mgiiams |?! | L24] 


and Pytte and Winsor For nitrogen, com- 
parisons will be made between the experimental results of 
this study and the Stine-Watson analytical solution. In 
making these comparisons, emphasis will be placed upon 
ascertaining the influence of nonequilibrium effects on the 
arc parameters. In comparing the anode heat transfer model 


(GRA 


of Pfender, et a with the experimental results of 





I2 


this study, attention is given to determining the ability 
of this model to accurately predict the anode heat flux, 


particularly in light of available empirical correlations. 


Iv.2 Comparisons For Argon 
Although the equilibrium model of Bower and Incro- 
pera and the nonequilibrium model of Clark have been ad- 
equately discussed previously, more must be said of the 
Stine-Watson analytical solution before uSing it in sub- 
sequent comparisons. As reported by Jedlicka and Seana 
the expressions for the electric field intensity and total 


wall heat flux which result from this closed form solution 


are given by 


o -CoL/tt 1/2 
E. aa [l-e ] (Ve) 
and 
F -CL/r 1/2 
q=C, > [l-e ] ; (IV.2) 
a 


where the constants Cor C3, and Cy are specific functions 
of the constant pressure gas transport properties, d 1s 
the duct diameter, and L is the duct inlet length. When 
applied to the asymptotic region (L + ~) of the constrict- 
ed tube arc, Equations IV.1 and IV.2 reduce to the expres- 


sions: 
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“3 
es ~ J (IV. 3) 
- il 
a= Cy 42 , (IV .4) 


Basing the above equations upon the tube radius and using 


the following results for the constants C3 and ous 


2.4 
C, = “ (IV.5) 
3 agi’? 
« 04363 ; 
Cy = agi? (IV.6) 


Equations IV.3 and IV.4 take the form 


mig 2 cua) 
and 
0.383 I 
q. = —= ° (IV .8) 
a Agi’? R- 


Ag (mho/watt) represents the assumed linearized relation- 
ship between the electrical conductivity and the heat flux 
potential ¢ (¢ = {kdT), and the units for the remaining 


variables are E. ,(volts/m), q, (watts/m*) , I(amp), and 
I 
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R(m). While the results reported by Stine and watson /+?! 


are primarily for nitrogen and air, this simplified model 


a7) 


has been modified by Bower to include the appropriate 


linearized property relationships for argon. By applying 
the value of Ag reported by Bower, the simplified Stine- 
Watson expressions (for argon with R = 0.5 cm) are then 


given by 


Ea = 9.55 volts/in (IV.9) 


and 
qo fey eeeeriin) = 63.4 I(amp) . (TV 31 0} 


The asymptotic values of the electric field inten- 


Sity (E ), total wall heat flux (qi), fraction of radia- 


Xx,a 


tive to total wall heat flux (g /4,) 1 mixed-mean en- 


rad,a 
Baalpy (h ), as computed 


[25] 


), and mean temperature (T 
[20] 


m,a m,a 


from the equilibrium and nonequilibrium models, 
are listed in Table IV.1 along with the corresponding ex- 
perimental results for an argon arc. The values from 
Table IV.1, along with results obtained from Equations 
IV.9 and IV.10 for the Stine-Watson model, were used to 
generate the plots Saown da Figures 27.1 through IV-.o. 


The large disparities (in both the magnitude and 


the functional dependence on arc current) between the 
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a Caragre IV .1) ae qs (Figure IV.2) distributions com- 
puted from the simplified Stine-Watson model and those ob- 
tained experimentally are thought to be the result of two 
oversimplifications inherent in the model. [In particular, 
the use of linearized property relationships and the as- 
sumption of negligible radiation losses, particularly at 
the higher arc currents, are thought to be the primary 
cause of these disparities. 

The comparison between the numerical solution by 
Bower and Incropera and the experimental results indicate 
that the equilibrium theory agrees favorably with experi- 
ment in the 100 amp range, while it overpredicts and under- 


predicts the experimental data for both E. and q, at 


,a 
lower and higher arc currents, respectively. A plausible 
explanation for the overprediction at the lower arc cur- 
rents can be obtained by considering the following expres- 
sion for Ohm's law 


ee (IV.11) 


x an f® orar 
O 
If the low current arc is characterized by significant de- 
partures from thermochemical equilibrium, the electrical 
conductivity used in the equilibrium solution will be less 


than the actual value of o in those regions where nonequi- 


librium conditions prevail. Consequently, the theoretical 
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Figure IV.1 Comparison Between Experiment and Theory 


for the Asymptotic Electric Field Inten- 
sity in Argon 
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evaluation of the integral of the electrical conductivity 
over the tube cross-section will be reduced, which, ac- 
cording to Equation IV.1l, will result in an overprediction 


et E increased Joule heating, and thus an accompanying 


x,a’ 
overprediciton of the total wall heat flux. While this 
nonequilibrium effect suitably explains the lower current 


comparisons, it does not explain why the equilibrium theory 


apparently underpredicts E.. a and q, at the higher arc 


currents. 
Recent theoretical studies by Kruger !?9! | Incropera 
and Piegase2>). and Clark !?>! and recent experimental 


[21] 


spectroscopic data reported by Bott and by Incropera 


and Pe riieacis '+-! 


indicate that arc nonequilibrium effects 
are most severe at low arc currents and diminish with in- 
creasing arc current. On this basis, the experimental re- 
sults and equilibrium predictions which are compared in 
Figures IV.1 and IV.2 should converge with increasing arc 
current when, in fact, they diverge for I > 100 amp. Ac- 
cordingly, the equilibrium underpredictions for Pe and 
ci at higher arc currents would not appear to be a result 
Of wall-induced nonequilibrium. 

A plausible explanation for the high current dis- 
parity is the inability of the equilibrium model to prop- 


erly account for conduction effects at the elevated temp- 


eratures associated with high current operation. In 
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particular, it 1s possible that the equilibrium model is 
underpredicting the energy transfer to the tube wall due 
to conduction. This contention is supported by the com- 
parisons for the fraction of radiative to total wall heat 
flux shown in Figure IV.3, as well as by the mean enthalpy 
and mean temperature comparisons of Table IV.1. Specif- 
ically, it 1s noted from Figure IV.3 that, while excellent 
agreement exists between experiment and the nonequilibrium 


predictions for gq /Q,) the equilibrium prediction, for 


rad,a 
I > 100 amp, is nearly double the nonequilibrium and ex- 
perimental distributions. In comparing the equilibrium 


and nonequilibrium predictions for q /A1 it may be 


rad,a 
argued that the observed differences are possibly a re- 
sult of a dissimilarity in the radiation source density 
terms used in the two solutions. This is not the case, 


[25] that, in the 


however, for it has been shown by Clark 
limit of equilibrium, the radiation model employed in the 
nonequilibrium solution agrees favorably with the radia- 
tion source density used by Bower and Incropera. More- 
Over, both of these plasma radiation models are consistent 
with the assumed wavelength distribution of the radiation 


incident upon the surfaces of the radiation gage (see the 


/q. is therefore 


Appendix). The overprediciton of q a 


rad,a 


most likely due to the higher temperature levels which 


Characterize the equilibrium solution. This fact is 
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further evidenced by the De comparisons of Table IV.1 as 


,a 
well as by the comparisons between equilibrium, nonequilib- 
rium, and experimentally determined temperature pro- 
files!) | 

On the basis of these arguments, it is concluded 
that, for the higher (I > 100 amp) arc currents in argon, 
the reduction of the equilibrium temperature profiles 
needed to bring about a corresponding reduction in the 
electrical conductivity and thus an increase in Ex a and 


qd (along with a decrease in g /Q, and ae a! could be 


rad,a 
achieved if the equilibrium model were to provide for 
greater energy loss by conduction. In an effort to con- 
firm this premise, some attention was given to the thermal 


conductivity model of pevoto!?"! 


which was employed in the 
equilibrium solution. 

In addition to the conventional translational mode 
of energy conduction, there will be a transport of chem- 
ical energy (in the form of the energy of recombination) 
through the ordinary (or concentration) diffusion of elec- 
tron-ion pairs towards the outer regions of the arc. 

While the global conservation equations employed in the 
equilibrium model do not explicitly account for ordinary 
diffusion, an attempt is made to implicitly account for 


the resulting transport of ionization (or recombination) 


energy through use of a reactive thermal conductivity 
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[372 


(A). From the work of DeVoto A, is defined by the 


R R 


heat flux vector expression 


Ge= aay Ta ARVT ; Funv eral) 


where k' is the translational thermal conductivity. This 
expression evolves from the more general form of the heat 
Flux vector equation by replacing the diffusion contribu- 
@ion to q by the term -A,VT, wherein it is assumed that 

the gradient of concentration can be replaced by a temper- 
ature gradient through the equilibrium law of mass action*. 
The point to be made here is that the critical assumption 
inherent in DeVoto's formulation of the reactive thermal 
conductivity is that of local chemical equilibrium. More 
specifically, it is assumed that, in the case of a chem- 
ically reacting mixture, all forward and reverse reaction 
rates are sufficiently fast to preclude any effect on the 
local composition due to ordinary diffusion. Otherwise, 
use of the equilibrium law of mass action to relate con- 
centration and temperature gradients could not be justified. 
Conversely, in a constricted tube plasma which is charac- 


8 
terized by strong concentration gradients, it is certainly 





* 
This derivation is treated in detail by Butler and 


[51] [52] 


Brokaw and by Fay 
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possible that the local atomic composition will be signif- 


Leeann 


Beantly influenced by ordinary diffusion 
case the assumption of local chemical equilibrium is in- 
valid. It is therefore suggested that the high current 
disparity between experimental and equilibrium theoretical 
results is due, in part, to a chemical nonequilibrium ef- 
fect on the reactive thermal conductivity. 

In connection with the above discussion, it is im- 
portant to note that this observation is seemingly incon- 
Sistent with previous considerations which indicated that 
the departure from thermochemical equilibrium is less se- 
vere at higher arc currents. However, in comparing DeVoto's 
calculations for the translational and reactive thermal 


[37] 


conductivities it is noted that the percentage contri- 


bution to the effective thermal conductivity (k = k'+h)) 


eer 


due to Ap increases from zero at T = 6,000°K to approxi- 
mately 50% at T = 14,000°K. Accordingly, any errors which 
may be introduced as a result of using DeVoto's formulation 
fOr Ap would become increasingly more significant at higher 
arc temperatures. In this respect, it is suggested that the 
failure of the equilibrium solution to provide for suffi- 
Cient energy loss by conduction at the higher are currents 

m= > 100 amp) is due to both its inability to allow for 


sufficient transport of ionization energy to the tube walls 


maa its inability to explicitly account for the effects of 
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diffusion on the local gas composition. In particular, it 
is suggested that the thermal conductivity values (specif- 
acally, Ap) employed in the equilibrium solution should 

be higher in order to bring the equilibrium theory predic- 
tions into better agreement with both the experimental 
observations and the nonequilibrium predictions. This con- 


[15] and 


tention is supported by the statements of Maecker 
by the comparisons reported by Morris, et al, (14! , These 
indicate a significant increase (by as much as a factor 
of two) in experimentally determined values of thermal 
conductivity over the corresponding predictions by DeVoto 
for temperatures greater than 10,000°K. 

While attempts have been made to interpret Figures 
IV.1 through IV.3 on the basis of nonequilibrium effects 
and thermal conductivity considerations, final judgement 
on the adequacy of the equilibrium model cannot be made 
until the combined experimental and theoretical uncertain- 
ities are considered. In an attempt to determine the 
range of uncertainties associated with the equilibrium 
solution, a parametric study was conducted !>3! wherein 


the properties (co, k and PR) were selectively varied 


EFF’ 


Over their respective limits of estimated accuracy 


iG, + 15%; en + 10%; and Pas + 15%)*. It was determined 


*The theoretical uncertainties for o and k were set by 


eee 


DeVoto (private communication with igeeeer while 


the uncertainty in Pr was arbitrarily sel@ceed >), 
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Beom this study that the upper uncertainty limit for E, - 
a 


and q, resulted from decreasing o by 15%, increasing KOPF 
by 10%, and increasing Pr by 15%, while the lower limit 


resulted from increasing o by 15%, decreasing k by 103%, 


Bae 
and decreasing Pr by 15%. The corresponding uncertainty 
bands for the equilibrium predictions are shown in Figures 
IvV.4 and IV.5 along with the experimental data and the 
related experimental uncertainities, Bad (+3%) and 

q, (46%). On the basis of the overlap between the theoret- 
ical and experimental uncertainty bands, and in view of 
the previous discussion concerning the use of DeVoto's 
formulation for the thermal conductivity, it is concluded 
that, to within the existing knowledge of plasma transport 
and radiation properties, single-fluid models are suitable 
for predicting asymptotic values for the electric field 
intensity and total wall heat flux, at least for the range 
of conditions in argon considered in this study. Converse- 
ly, the equilibrium model is not suitable for predicting 

ad ,a/ ta over the entire range of conditions investigated, 


and T at the 
m,a 


a 


nor is it adequate for predicting ha Z 


a 


lower (< 50 amp) and higher (> 150 amp) arc current con- 
ditions. In addition, it is felt that the equilibrium 


model would not be suitable for predicting EY and qo 


7a 


if the above comparisons were extended to include condi- 


tions of lower and higher arc currents. In this respect, 
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Figure IV.4 Comparison Between the Combined Experimen- 


tal and Theoretical Uncertainties for the 
Asymptotic Electric Field Intensity in 
Argon 
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it is concluded that a more rigorous theoretical treat- 


ment, such as the nonequilibrium solution. by clark !27! , 


would be required to accurately model the asymptotic re- 
gion of the constricted tube argon arc. 

Before concluding this section on comparisons for 
the asymptotic region in argon, some explanation is in 
order concerning the lack of agreement between the exper- 


imental value and the nonequilibrium prediction for h., a 


at I = 35 amp (Table IV.1). Calculation of ho = in the 


oat 


nonequilibrium solution is effected through the expres- 


Bion |>4! 


; fe opuhrdr 


h (iV273) 


UG fe ourdr 


where o(r) is the local density of the plasma and u(r) is 
the global axial velocity component. The radial enthalpy 
distribution, h(r), required for the above calculation is 
given by 

Kae ion 


[aT (r) + T, (r)J] + 
My = h My 


h(r) = Coie ee, (IV.14) 





rj 


where the first two terms represent the contributions due 
to the electrons and heavy particles, respectively, and 
the remaining term represents the chemical enthalpy of 


recombination (the specific variables are defined in the 





oul 


nomenclature). In formulating Equation IV.14, the addi- 
tional enthalpy contributions due to internal electronic 
excitation and radiation source density effects were as- 
sumed negligible. The governing equations for the asymp- 
totic region are obtained from the general set of conser- 
vation equations by setting all axial derivatives (with 
the exception of the pressure gradient, dp/dx) equal to 
zero !25! | This manipulation, along with the assumptions 
of negligible viscous dissipation and pressure work, re- 
sults in a decoupling of the hydrodynamic and thermal flow 
fields. Consequently, the velocity profiles required in 
the above definition for hiss _ can only be obtained by 


f 


solving the following asymptotic axial momentum equation 


Q) 


du, _ dp _ 
gy (ur gel - a = 0 - ee 


KI 


It is the uncertainty associated with the determi- 
nation of dp/dx in Equation IV.15 that renders the non- 
equilibrium solution unable to accurately predict u(r) 
and thus Dna at arc currents less than 50 amp. For the 
same reason, nonequilibrium predictions for ea at 
currents greater than 100 amp are unrealistic - the agree- 


ment in the 50 to 100 amp range is thought to be fortu- 


ZEOUS. 
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IV.3 Comparisons for Helium and Nitrogen 


Asymptotic values for helium of the electric field 
intensity, total wall heat flux, mean enthalpy, and mean 
temperature computed from the equilibrium model by Bower |>>} 
are listed in Table IV.2. Also listed in this table are 
the experimental data for helium reported by McKee, et 


mp, 132] oie as well as the 50 ampere data of 


and Emmons 
this study. The data from Table IV.2 were used for the 
plots of Figures IV.6 and IV.7, wherein it is noted that 
the experimental qd values (with the exception of the 
single qo value measured calorimetrically in this study) 
were obtained from the asymptotic energy balance equation 
(Equation III.2) and the experimental Bes data. 

In comparing the equilibrium theory predictions 
With the experimental results, it is readily apparent 
that, unlike the previous comparisons for argon, the 
Single-fluid model is totally deficient in its attempt 
to accurately describe the overall energetics of the heli- 
um arc. The reason for this deficiency is believed to be 
a result of the severe nonequilibrium conditions which 
are known to exist in low current (I < 400 amp), atmo- 
spheric helium arcs. 

The pronounced departure of helium arc plasmas 
from conditions of local thermodynamic eqiulibrium has 


been observed and discussed by others, among the first of 
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[12,13,50,56] 7hese 


which were investigators at Harvard 
studies were primarily concerned with determination of 
the plasma transport properties, o and k, through use of 
the so-called inversion method. As discussed by Emmons 


and anaes) 


this method, when formulated under the as- 
sumptions of local thermodynamic egqiulibrium and negli- 
gible radiative losses, dictates that the reduced arc- 
column characteristics (EY al versus Ey ak should be 
related by a universal curve independent of tube radius. 
However, in performing electric field intensity measure- 
ments on helium arcs of varying tube radius, Emmons !13] 
determined that the parameters Exar and Ey ak were not 
uniquely related. Since the assumption of negligible 
radiation is apparently valid for low current (I < 400 


[57] Emmons concluded that 


amp), atmospheric helium arcs 

helium arcs must be characterized by a significant de- 

parture from conditions of local thermodynamic equilibrium. 
Further evidence of severe nonequilibrium conditions 


[50] Bott (6! 


in helium arcs is provided by Kruger , and 
Aleksandrov, et ai, '>8! | Specifically, these investigators 
concluded that the free electrons and upper excited 

atomic states, while in equilibrium with one another, are 
not in equilibrium with the lower excited (n < 3) and 
ground states. In addition, Kruger notes that, in the 


outer regions of the arc, the electron number densities 


95 


may be several orders of magnitude higher than would be 
predicted if the local composition were determined from 
the assumption of Saha equilibrium. The existence of 
Significant departures from equilibrium can also be in- 
ferred from the helium data acquired in this study. As 
indicated in Table IV.2, energy balance measurements per- 
formed on the helium arc dictate an asymptotic mean tem- 
perature of 1350°R (or 750°K) and yet, even at such a 


, there is sufficient ion- 


comparatively low value of Ta x 
a 


ization to produce a luminescent arc which is capable of 
Sustaining a current of 50 amperes. This observation 
Suggests the existence of a severe thermochemical non- 
equilibrium condition in the helium arc. Similar ob- 


[58] whose 


servations were made by Aleksandrov, et al. 
spectroscopically determined electron temperatures were 
nearly 1.5 orders of magnitude greater than their mean 
temperatures based upon energy balance measurements. 

On the basis of these findings, the comparisons 
between experiment and the equilibrium theory of Bowers 
(Figures IV.6 and IV.7) are virtually meaningless. That 
is, since the plasma tube flow situation which the single- 
Fluid model is attempting to describe bears little re- 


semblance to the severe nonequilibrium state which actu- 


ally exists, a systematic discussion of the differences 
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is not possible. As indicated in Figure IV.6, the pre- 
@eetion of the arc-column characteristics in helium by 


Pytte and Williams /48! 


offers some improvement over the 
numerical solution by Bower; however, this theoretical 


model continues to underpredict E by some 20%. This 


X,a 
underprediction by the Pytte-Williams model is most likely 
due to the use of equilibrium calculations for o and «1391 
as well as to use of the Elenbaas-Heller equation. Pytte 


[24] 


and Winsor have given greater consideration to non- 
equilibrium in atmospheric helium arcs by proposing a two- 
temperature theory which is based upon the asymptotic 
moment equations of continuity, momentum, and energy as 
derived from an electron transport (Boltzmann) equation. 
While this model accounts for thermal nonequilibrium be- 
tween the electrons: and heavy particles, it suffers from 


ae and 


use of equilibrium thermal conductivity values 
use of the Saha equation (evaluated at the local electron 
temperature) to determine the local composition. In ob- 
serving the excellent agreement which exists between the 
Pytte-Winsor predictions and Emmons' experimental data for 
Ba (Figure IV.6), it would appear that use of the Saha 
equation and egqiulibrium thermal conductivities is justi- 
fied. This agreement is misleading, however, for Pytte and 


Winsor adjusted their electrical conductivity, 9, in order 


to obtain the best fit to Emmons' experimental reduced 





wo 


arc-column data for p = 1.0 atm and R = 0.5 cm. [In apply- 
ing the resulting value of o to calculations for tube 
radii other than R = 0.5 cm, the Pytte-Winsor model does 
predict the trends indicated by the experimentally de- 


termined shifts in the E, ah versus ES al curves; however, 


ga a) Ee 


it fails to predict the magnitude of these shifts 
is believed that this deficiency results from the failure 
of the Pytte-Winsor model to properly account for chemical 
nonequilibrium. Accordingly, the need remains for a rig- 
orous, theoretical, multifluid model which treats the 
thermochemical nonequilibrium flow of helium in the inlet 
and asymptotic regions of the constricted tube arc. 

San lar em, the: author knows of no multifluid 
analyses in the literature from which predicitons of the 
global parameters, Ey and q can be made for the inlet 
and asymptotic regions of nitrogen arcs. Consequently, 
the following comparisons for nitrogen are limited to the 
Single-fluid analytical model by Stine and watson (18) . 
As was the case for argon, the generalized Stine-Watson 
expressions for ae and qd. (Equations IV.7 and IV.8) 
must be adapted to the particular gas and tube radius 
being considered. For this purpose, the values of E 


[3] 


and R for nitrogen reported by Shepard and Watson were 


xX,a 


used with Equation IV.7 to determine Ag. This value was 
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then used to give the following simplified expressions for 
the Stine-Watson model, as applied to the asymptotic re- 


gion of a nitrogen are with R = 0.5 cm. 


Era = 26.4 volts/in (IV.16) 


i. = 7 Ses Sis ak (LV) 


These relationships, along with the experimental 
data from Figures III.6 and III.9, were used for the com- 
parisons of Figures IV.8 and IV.9. Relative to previous 
comparisons for argon, application of the Stine-Watson 
model to nitrogen provides better agreement with exper- 
iment. While the temperature dependence of the radiation 
source densities (P,) for argon and nitrogen are nearly 


the same !°9] | 


the mean temperatures in nitrogen are some 
590% lower than the aomnemaenelsins values of ee Le angon 
(Figure III.13). Therefore, the percentage contribution 
to the total wall heat flux due to radiation from a nitro- 
gen arc plasma would be considerably less than the corre- 
sponding contribution for an argon plasma of the same arc 
current. As such, the negligible radiation assumption 
implicit in the Stine-Watson model is less severe for 
nitrogen; nevertheless, as shown in Figures IV.8 and IV.9, 


this assumption continues to render this model inaccurate 


at higher arc currents. 
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In making further comparisons between the three 
gases, it is noted from Figure III.6 that, for a given arc 
current, successively higher voltages are required to 
sustain arc operation in going from argon to nitrogen to 
helium plasmas. Moreover, the increased values of EY a! 
which reflect the dissociation energy of No and the rela- 
tively high ionization potential of He, necessarily result 
in the proportionately higher values of qo for nitrogen 
and helium which are indicated in Figure III.9. On the 
basis of the increased values of total wall heat flux for 
nitrogen and helium, it might be expected that the mean 
temperatures of these two gases would be higher than the 
corresponding values of a for argon. However, this is 
Clearly not the case. As indicated in Figure III.1l3, 
the Dy es values for nitrogen and helium are considerably 


g 


lower than the Tis 7 values for argon. In fact, at I = 50 
g 


amp, the experimentally determined value of oe fOr 
helium is more than an order of magnitude lower than the 
corresponding value for argon, and yet the calorimetri- 
cally measured value of qa for He exceeds the correspond- 
ing value for Ar by a factor of nearly six. While this 
behavior is seemingly anomalous to the rules of convection 
heat transfer, it is recognized that any attempt at ex- 


planation on the basis of the global measurements of this 


study would be fatuous. While a partial explanation for 
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the increased wall heat transfer for the helium arc over 
that for argon can be inferred from the characteristically 
higher values of thermal conductivity for He (e.g., at 


[37,40], 


m= 5,000°K, Kure = a complete explanation 


9 Kar 
cannot be gained without consideration of the microscopic 
transport phenomena which characterize these gases. 

The order-of-magnitude difference between the mean 
temperatures in argon and helium can be explained on the 
basis of the relative specific heats for these two gases. 
That is, while the experimentally determined mean en- 
thalpies for He and Ar are nearly equal (Tables IV.1 and 
IV.2), the specific heat for helium is nearly an order of 
magnitude higher than the corresponding value for argon. 
When evaluated at the experimentally determined values 
|f Dna for argon and helium (Tables IV.1 and IV.2), the 
ratio of the constant pressure specific heats (He to Ar) 


Me 9.113839] 


IV.4 Comparisons for the Anode Region 


Before comparing experimental and theoretical 
anode heat transfer results, it is necessary to discuss 
the theoretical model and its relevance to measured 
quantities. Since the major portion of the experimental 
data acquired in this study is for argon, the following 


discussion and comparisons will deal only with this gas. 
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An anode heat transfer model which is frequently 


found in the literature (6"/ 126-30) is one based upon an 


energy balance over the surface of the anode. 





Anode Heat Balance Model 


Figure IV.10 


With reference to Figure IV.10, this balance is of the 


form 


; (IV.18) 
rt+c e 


where 


Qn = rate of energy loss to the anode cooling 


water, 





"Or 
tl 


ier Cc 


Oe 
tt 


Q5b1 


On the 
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rate of energy transfer to the anode sur- 
face by convection and radiation from the 


plasma column, 


rate of energy transfer to the anode sur- 


face by electrons, 


rate of energy transfer from the anode sur- 
face due to ablation of the anode material, 


and 


radiation emissive power of the anode sur- 


face. 


basis of observed negligible degradation of 


the anode surface after some 60 hours of operation and 


from calculations for a black surface emitting at the 


melting temperature of copper, the last two terms in Equa- 


tion IV.18 are known to be insignificant compared to 


calorimetrically measured values of on and can therefore 


be neglected. 


In addition, the energy transferred to 


the anode surface by the electrons is commonly obtained 


from the expression 


meme 


[6] 


2 


: as kT. + U,) J oni, (IV.19) 
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where 

io.5 = 

= 3 kT. = thermal energy of the electrons as they 
enter the region of "anode fall", 

IU, = additional energy acquired by the elec- 
trons due to acceleration through the 
anode fall region, 

Io. = "condensation" energy of the electrons 
penetrating the anode surface (>. = 
work function of the anode material), and 

E = fraction of energy retained by the elec- 


trons after collisions with heavy parti- 


cles within the anode fall region. 


In addition to assuming negligible OLan and Our 
Pfender, et ay lGr71 in working with Equations IV.18 and 
IV.19 assume a value of unity for € and assume that Un 
may be represented by the difference between the anode 
potential (measured with respect to the cathode) and the 
potential obtained by linearly extrapolating the fully- 
developed potential gradient (EL a! to the location of 
the anode (see Figure IV.11). In using these assumptions 


with selected values of T, = 10,000°K and on = 3.7 volts, 


Equations IV.18 and IV.19 reduce to 
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Q, = Q 


+ I (5.85 + Uy) - (IV.20) 


eC 


If the anode is preceded by a thermally fully-developed 
region and if the mean enthalpy of the plasma at the anode 


remains close to h, ., Pfender, et al.argue that ore 


tf 


becomes the asymptotic value of the total wall heat flux 


(q.)- Accordingly, Equation IV.20 takes the form 


O + 26.6 I (5.85 + Un), (iV ee) 


A 4a 


where the constant 26.6 accounts for the area of the anode 
surface (Qn and qs have units of ipo ee oD and the 
appropriate conversion factor (0.0569 oa: 

Using the experimental data of this study for qa 


and U a comparison is made between the values of Qn 


a! 
calculated from Equation IV.21 and the measured values of 
Enis study (Figure III.15). The results of this comparison 
are shown in Table IV.3, wherein it is noted that the 
theoretical model overpredicts the corresponding experi- 
mental values by a factor of approximately two. It is felt 
that a portion of this rather gross overprediction is due 
to errors associated with assuming that € = 1.0 and 


assuming that U, can be determined by the aforementioned 


A 
linear extrapolation procedure. These errors are made 
evident by the anode energy balance measurements of this 


study which reveal a significant increase (~ 30% above 


0:9 
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the value of Ra in the enthalpy of the gas as it passes 
through the anode region. In order to obtain such an 
increase in enthalpy, there must occur a significant trans- 
fer of electrical energy to the gas in excess of that 
required to maintain the asymptotic conditions ee and 
thus q.)- The electrons, therefore, will in fact lose a 
portion of their energy (IU,) to the heavy particles in 

the fall region and, consequently, & will be considerably 
less than unity. The failure of the theoretical model 
(Equation IV.21), when applied to the anode configuration 
of the constricted tube arc, to account for this additional 
transfer of electrical energy 1S perhaps best understood 
by considering a typical anode region potential distri- 
bution. With reference to Figure IV.11, the theoretical 
model assumes that, of the total anode potential rise, En? 
the portion, FE, Uy is required to maintain the asymptotic 


value of ne e and the remaining portion, U results in 
? 


na 
energy transfer, via the electrons, to the anode surface 


(l.e., IU,). In contrast, the experimental results 


indicate that, in addition to that portion of the total 


anode potential rise required to maintain h a portion 


m,a’ 


of EA ts aso Seeded" te sustain the increase in the 


enthalpy of the gas, Ah as it passes through the anode 


ie 


region. Consequently, a fraction of the energy acquired 


ee 


— 








Figure IV.1l Partitioning of Total Anode Fall 
Potential 


by the electrons (IU, ) is in fact transferred to the gas 
through electron-heavy particle collisions and cannot 
contribute directly to the total anode heat transfer. 
Hence, it is concluded that the overpredictions for Qa 
which result from the anode heat balance model are at 
least partially due to the combination of assuming that € 
is unity and that the anode fall potential is given by 
the linear extrapolation procedure. 

To eee consider these assumptions as the major 


sources of error in Equation IV.21, the data of this study 


pa? 


were used to determine both the potential drop required 
to sustain the gas enthalpy rise Ah, and the potential 
drop, to be designated Un , which results in the transfer 


1). 


of electron energy directly to the anode surface (IU y 


In this connection and with reference to Figure IV.1l, it 


is assumed that the total anode potential rise, E can be 


a 
subdivided into the fraction of this potential required 
to sustain the asymptotic conditions, E7U,/E, i the frac- 
tion of this potential required to sustain the anode re- 
gion enthalpy rise (Ah,), U,- UE, ; and the fraction of 
this potential resulting in energy transfer to the anode 
surface by the electrons, Un /E,- According to this sub- 


division, and retaining the values of T. and dn used in 


Equation IV.21, energy balance considerations require that 
Qn = q, = 2066 -beeoe co oa Un) ; (IV .22) 
ema that 


I (U,. - U_,) = mAh, , (i722 3) 


A! 


where, again, Ah, is the enthalpy acquired by the gas in 
traversing the anode region. The anode heat flux predic- 
tions which result from calculations performed with this 
modified theoretical model are presented in Table IV.4 


along with the experimentally determined values of Ah, - 
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While the values of U computed from Equation IV.23 


A 
are significantly lower than those given by the linear 
extrapolation method (Un), their use in Equation IV.22, 
along with experimental values of qa! still result in 
sizeable overpredictions for the total anode heat flux. 
It appears then that more appropriate values for T. and 
on must be determined in order to bring the theoretical 
predictions into better agreement with experiment. These 
additional considerations would necessarily require further 
knowledge of electron temperature distributions in the 
Vicinity of the anode as well as a more rigorous consid- 
eration of the plasma-electrode surface interaction. Ac- 
cordingly, the use of this analytical model to predict the 
total heat transfer to the anode of a constricted tube 

arc appears to be of questionable value. This is espec- 
ially true when the model is compared with existing em- 
pirical correlations, such as those appearing in Figure 
III.15. In connection with these empirical correlations, 
it is noted that in cases where the anode is not preceded 
by a fully-developed region, it is necessary to include 

a dependency on arc length and operating gas flow rate, 
[29] 


as is the case, for example, of the Eberhart and Seban 


anode heat transfer correlation. 
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CHAPTER V SUMMARY AND CONCLUSIONS 


This work provides experimental results for the 
axial distribution of the electric field intensity and 
total wall heat flux for the inlet, asymptotic, and anode 
regions of atmospheric argon, helium, and nitrogen arc 
plasmas. Whereas data are presented for a wide range of 
laminar flow rates (0.03-0.11 l1b/min) and are currents 
(35-240 amp) in argon, the amount of data presented for 
helium and nitrogen arcs is considerably less extensive. 
The range of operation for these two gases was restricted 
by the power supply employed in this study. In addition 
to the development of the diagnostics required for the 
electric field intensity and total wall heat flux measure- 
ments, considerable effort was given to the development 
of a unique gage to measure the radiative component of 
the total wall heat flux. This instrument has been used 
successfully in providing the radiative wall heat flux 
data which are presented for the asymptotic region in 
argon. The accuracy of the experimental measurements, 
as well as the consistency of the arc facility, is assessed 
on the basis of the observed precision of local energy 


balance checks. Further confidence in the experimental 
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measurements is inferred from the favorable agreement 
which exists between the data of this study and the experi- 
mental results reported by other investigators. 

Through comparisons between the experimental data 
and available equilibrium and nonequilibrium numerical 
solutions, an attempt is made to ascertain the adequacy of 
single-fluid models to accurately describe overall arc 
energetics. On the basis of the measurements themselves 
and the comparison with theoretical results, the following 
conclusions were derived. 

l. The increase in duct wall potential with current 
drawn from the segment is primarily a "cold" 
boundary layer effect and, therefore, the 
"floating" duct wall potential provides a more 
accurate measure of the actual plasma potential, 
at least for the copper wall-argon plasma 
combination employed in this study. 

2. Variations in the duct entrance conditions can 
have a significant effect on the experimentally 
determined and fonecic cnn predicted electric 
field intensity and total wall heat flux dis- 
tributions for the inlet region of the con- 
strictor. As a result, meaningful comparisons 
between experiment and theory can be made only 


for the asymptotic region. 
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Closed form analytical solutions which are 

based upon the use of linearized property 
relationships and which assume negligible 
radiation are inadequate for predicting electric 
field intensity and total wall heat flux. This 
1s particularly evident at the higher current 
levels in argon and nitrogen. 

TO within the combined experimental and theoreti- 
cal uncertainties, the numerical solution of a 
rigorous single-fluid model is suitable for use 
in predicting the electric field intensity and 
total wall heat flux for the asymptotic region 
in argon, at least over the range of conditions 
investigated in this study. The equilibrium 
solution does not, however, appear suitable 

for predicting the asymptotic mean temperature 
or the radiative wall heat flux over the same 
range of conditions in argon. There is a small, 
but discernable, nonequilibrium effect at the 
lower arc currents in argon, as a result of 
which the equilibrium solution tends to over- 
predict both the electric field intensity and 
total wall heat flux for currents less than 

100 amp. Moreover, the tendency of the equili- 
brium solution to underpredict these two para- 


meters at arc currents in argon greater than 
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100 amp is believed to be due to a chemical 
nonequilibrium effect which renders the 
equilibrium solution deficient in providing 
for sufficient energy loss to the tube walls 
by conduction. 

5. Solutions which properly account for thermo- 
chemical nonequilibrium, as well as for all 
pertinent transport phenomena, provide pre- 
dictions for the electric field intensity and 
the convective and radiative wall heat transfer 
which are in excellent agreement with the 
experimental results over the entire range of 
asymptotic conditions investigated in argon. 

6k Atmospheric helium arcs are characterized by a 
Significant departure from local thermodynamic 
equilibrium (at least for arc currents less 
than 250 amp) and, as a result of these 
departures, single-fluid models are totally 
inadequate in representing their thermal state. 

With reference to the last three of the above con- 

clusions, it is presumed that further comparisons between 
experiment and equilibrium and nonequilibrium solutions, 
if extended to include are conditions characterized by 


higher order nonequilibrium effects, would lead to an 
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improved understanding of nonequilibrium arc physics. 
Accordingly, it is recommended that the present arc 
facility be extended to permit stable, low current 

(I < 35 amp) operation in argon and to permit stable 
operation over a greater range of operating conditions 
in helium and nitrogen. It is also recommended that the 


[20] 


solution by Bower and Incropera be expanded to provide 


equilibrium predictions for nitrogen arcs and that the 


solution by Clark 25! 


be appropriately extended to afford 
nonequilibrium arc predictions for helium and nitrogen. 

In addition to the above conclusions for the inlet 
and asymptotic regions, further comparisons between exper- 
iment and theory for the anode region indicate that the 
large (50-100%) overpredictions for the total anode heat 
flux which result from use of a heat balance model are 
primarily due to the relatively large uncertainties in- 
volved in specifying associated values of the physical 
parameters (Uns Tor and on) > Consequently, the theoreti- 
cal model is thought to be ill-suited for providing de- 
Sign estimates for the total heat flux to the anode of 
the arc constrictor, particularly in comparison with 
estimates which are based upon empirical correlations. In 


this respect, further experimental and theoretical study 


of plasma-electrode surface interactions is recommended. 
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APPENDIX RADIATION GAGE 


A.1 General Remarks and Review of Experimental 
Studies Relating to Radiative Heat Flux Measurements 


Since thermal radiation is often a primary mode 
of energy transfer in a constricted arc plasma, knowledge 
of the radiative contribution to the total wall heat flux 
is of concern in both the engineering application and 
the laboratory use of plasma arc generators. In addition, 
Beeowledge of the radiative wall heat flux is of value in 
obtaining an improved understanding of arc physics and 
plasma flow modeling. Accordingly, a portion of this 
study was directed towards the design, fabrication, and 
use of appropriate instrumentation to measure the radia- 
tive component of the total heat flux to the wall of an 
mec CONStrictor. 

While there are numerous studies in the literature 
dealing with the measurement of plasma radiative proper- 
ties, the author knows of no experimental attempt to 
deal explicitly with the measurement of the radiative 
and convective components of the total wall heat flux for 
an arc constrictor. There are described in the litera- 


ture, however, two diagnostic techniques which would appear 
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applicable to this particular measurement. The first of 
these involves a total radiation detector which utilizes 

a thin-film resistance sensor. This sensor consists of 

a thin metallic strip (usually platinum) which is depo- 
Sited on a substrate material (pyrex, quartz, etc.) and 
has been used primarily for transient, short duration 
measurements. The metallic strip acts as a passive re- 
Sistance element in an electrical circuit designed to 
record the variation of current through the sensor upon 
its exposure to a step input of radiation. From know- 
ledge of the temperature coefficient of resistance of the 
metallic element and the thermal conductivity of the sub- 
strate, the heat flux to the substrate can be inferred 

by analyzing the time variation in the surface temperature 
aS indicated by the electrical output from the sensor. 
Further discussion of this diagnostic technique is provided 


[61] 


by Reece and Theoclitus A major advantage of this 


method is that the calibration and data reduction procedures 
: . hoz On) 
are relatively straight forward and well established 4 
While the time limit of sensitivity (saturation 
time) for the thin-film sensor is typically on the order 
of 100 msec, this technique can be modified for use with 
steady state radiation sources. As discussed by Hardy 


[64] 


and Paddock , this feature is made possible by placing 
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a fast-acting shutter mechanism between the radiation 
source and the sensor. This modification would, of course, 
be required if the thin-film detector were to be used for 
the radiation measurements of this study. However, here- 
in lies the major disadvantage of this method. In parti- 
cular, the imposition of interface components (shutter 
mechanism, etc.) would place the sensor sufficiently far 
from the plasma source to require detailed consideration 
of the spectral sensitivity of the resulting cavity walls. 
A second device has been employed by researchers 
at Harvard to measure the total radiation from a steady 
state, constricted tube arc.* This device uses a re- 
cessed sensing element, across which a temperature 
difference is induced which is proportional to the ab- 
sorbed radiation. The sensing element consists of a 
blackened constantan shim which is positioned in a cir- 
cumferential (with respect to the axis of the cylindrical 
duct) slot. The slot walls are silver coated in order to 
reflect a greater percentage of radiation to the sensor, 
and the sensing element is recessed 1.5 inches in an 
attempt to minimize convective contributions to the 
measured heat flux. Calibration of the radiation meter 


was attempted by comparing the detector output against 


* A comparison is made in Chapter III between the Harvard 
results and the radiation data obtained in this study. 
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several radiation standards, which include filament and 
gas-filled lamps and globar units. 

This diagnostic technique is subject to several 
sources of error for which compensation cannot be entirely 
made by the calibration procedure. Consequently, accurate 
use of this device is contingent upon the successful 
isolation and evaluation of these sources of error and the 
subsequent introduction of appropriate correction factors. 
This matter deserves further discussion, for it serves to 
illustrate the problems associated with radiative heat 
flux measurements, particularly with techniques which 
employ recessed sensors. 

As discussed by aout certain correction 
factors must be applied to the experimental results ob- 
tained by use of the aforementioned device. This is 
necessary in order to account for the geometrical differ- 
ences between the conditions of calibration and actual 
measurement. Specifically, account must be taken of the 
difference between the relative positions (with respect 
to the sensor location) of the radiation standard and the 
arc. In addition the variation in geometry which results 
from the fact that the calibration standard (globar unit) 
is a surface emitter, whereas the arc is a volume emitter, 


must be considered. These correction factors may be 


ales? 1. 


obtained from geometrical considerations relating to the 
Calibration and actual measurement configurations. How- 
ever, aside from the geometrical effects, there remains 
a possible source of error due to the spectral sensitivity 
of the component parts, particularly the slot walls. If 
the reflectivity of the slot walls varies appreciably over 
the wavelength range of the incident plasma radiation, and 
if this variation is to be incorporated within the cali- 
bration of the detector, then the radiation standard must, 
of course, have the same spectral intensity distribution 
as the plasma. Unfortunately, such calibrated sources 
are not available. 

Another deficiency inherent in this radiation gage 
is the fact that the higher plasma radiation intensities 
could only be determined from a linear extropolation of 


[47) mat is, the detector output 


the calibration data 
during actual measurement was a factor of approximately 
six higher than the maximum calibrated detector output 

(as limited by the burnout temperature, 1625°C, of the 
globar unit). Techniques employing recessed sensors are 
therefore calibration limited in the sense that if the 
possible variations in the spectral response of the cavity 
walls (as well as the full output range of the detector) 


are to be incorporated within the detector calibration, 


the calibration source must have the same radiation 
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intensity, in both magnitude and spectral distribution, 


as the plasma for which measurements are desired. 


A.2 Radiative Heat Transfer Gage - General 
Design and Operation Considerations 


From evaluation of the aforementioned measurement 
techniques, three primary constraints were initially 
imposed upon the proposed detector design. First, the 
detecting unit must not alter the geometry of the con- 
strictor and, in particular, the use of recessed sensing 
elements should be avoided. Secondly, the detector must 
not be calibration limited in the sense that generation 
of its output characteristic is dependent upon an initial 
calibration against a standard radiation source. Thirdly, 
the presence of the radiation gage should not significantly 
alter local conditions in the plasma. This is necessary 
if useful comparisons between theory and measurements 
are to be realized. 

To adhere to the above criteria, it seemed rea- 
sonable to attempt to alter the wall surface characteris- 
tics to permit isolation and measurement of the radiative 
heat flux. The use of a transparent "window" affords a 
convenient means of isolation; however presently avail- 
able optical materials cannot provide uniform trans- 
mission over the entire spectral range of possible plasma 


radiation. In addition, the window would have to be 


L3$ 


thermally isolated from the sensor so as not to contri- 
bute to the measured heat flux through conduction. 
Further insight into a possible technique for de- 
lineating between the radiative and convective components 
to the total wall heat flux was gained from consideration 
of the optical characteristics of the duct wall surface. 
Expressing the radiative component in terms of the mono- 
chromatic absorptance (a, ) and the plasma irradiation 


(H)), the total wall heat flux takes the form 


So dconv i: drad 


oO 
pee | Ronee (A.1) 
O 
—— 4, . 
The emissive power of the duct wall surface (eoT, ) is 
neglected, since its magnitude is small in comparison to 
the other terms of the equation. Assuming, for the moment, 


that a, is known and that q may be determined by calori- 


r 

metric measurement, there remain two unknowns in Equation 
: unknowns 

ae lig, WE oysns7 and Hy The appearance of these two 

suggests the use of two surfaces (which differ in spectral 


sensitivity), in which case the application of Equation 


A.1 for surface-l yields 
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= + 
4) doonv,1 Joy iy al a (A. 2) 
and, for surface-2, 
co 
42 ~ Aconv,2 4 Pf ee Ss oo 


If both surfaces are located at the same axial location 
in the duct, it may be assumed that both the incident 
irradiation and the convective heat transfer to each sur- 
. = = H a 
face are the same (i.e., Hy 4 a, \ nd q 


= ° U , . 
Iconv ,2 Seon?) nder these assumptions, Equations 


conv,l 


A.2 and A.3 can be equated on the basis of equal convection 


to give 


co 
ae aaciey Jo H. da - Jay 2 H, da. - 


Ne a 


To achieve the largest possible difference between 
the calorimetrically measured values of the total wall 
heat flux (a, - Go) and thus to realize the maximum 
experimental resolution, the difference between the radia- 
tion absorption terms in Equation A.4 should be maximized. 
This suggests the use of a highly absorbing (black) sur- 
face in combination with a highly reflecting surface. 

For the optically black surface, say surface-l, it is 


assumed that suitable coating materials are available for 


which the monochromatic absorptance is independent of 
wavelength over the entire spectral range of the plasma 
irradiation. For such a surface, the radiative wall 


heat flux component reduces to 


J Oy H, d\ = Oy JH) di = a4 H (A 


Similarly for the reflecting surface (surface-2), 


r 2 


O a 


With these additional assumptions, Equation A.4 reduces 


to 


fo, , H, dA =a, fH) dk = aH . (A. 
0 


13 


o>) 


a) 


Knowledge of the absorptance values, rT and Qos along with 


the experimentally determined values of qy and qo will 


permit evaluation of the plasma irradiation. This value 


of H can then be employed in either of the following 


expressions to yield the value of the convective wall heat 


transfer: 


= + H A. 
q,-4 Oy ( 


or, 


He. (A. 


8 ) 


2) 
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To determine the radiative wall heat flux to the standard 
(copper) duct segments, this value of cee 1s subtracted 
from the total wall heat flux results (obtained at the 
same Operating conditions and axial location) determined 
from previous calorimetric measurements. 

As shown in Figure A.1, the radiation gage consists 
of two semicylindrical copper segments which are thermally 
isolated from one another, individually water cooled, 
and supplied with individual AT-transducers. An optical 
black (Nextel 101-C10 by Minnesota Mining and Manufactur- 
ing Company) coating is applied (and then baked) to the 
interior cylindrical surface of one of the half-segments, 
while the opposing half-segment is coated with a highly 
reflecting, thermally evaporated aluminum film.* The 


overall dimensions of the assembled radiation gage are 


precisely the same as those of a standard duct segment. 


A.3 Discussion of the Parameters and 


Assumptions Affecting the Operation 
of the Radiation Gage 


While the above radiation gage design appears to 
Satisfy all of the initially posed constraints and 
specifications, the following two basic assumptions re- 


quire further discussion: 


* The fabrication and optical properties of the reflecting 
surface are discussed in Section A.4. 
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1. The presence of the radiation gage does not 
Significantly alter the thermal state of 


the plasma in the vicinity of the gage and, 


2. The reflecting surface is characterized by a 


constant value of absorptance, a over the 


Den 
entire spectral range of the incident plasma 


irradiation (Hy). 


The validity of the first assumption is critical 
to the success of this diagnostic technique. That is, 
if the pressence of the highly absorbing and reflecting 
surfaces results in a significant change in the local 
thermal state of the plasma, then the assumptions of 
equal convection and irradiation required by Equation 
A./7 would be rendered questionable. Moreover, any such 
change in the local thermal plasma state would preclude 
combining the convective wall heat flux (determined from 
use of the radiation gage) with the total wall heat flux 
data (obtained from calorimetric measurements employing 
the standard, copper-wall duct segments) to determine 
the radiation flux. 

Wall heat transfer measurements were made at the 
radiation gage segment and the eight adjoining standard 


segments (four upstream and four downstream) in an attempt 





uy 


to experimentally determine if there existed any effect 

on the plasma radiosity distribution that could be 
attributed to the spectral characteristics of the gage 
surfaces. This experimental investigation in the 
asymptotic region of an argon arc was performed using both 
black-versus-reflector and reflector-versus-reflector 
surface combinations. The results showed no discernable 
difference between the total wall heat flux to the seg- 
ments adjoining the radiation gage and the qs values which 
were measured in the absence of this gage. Moreover, the 
total wall heat flux to both of the reflecting surfaces 
for the case of the reflector-versus-reflector combina- 
tion compares favorably with the measured total wall 

heat flux to the reflecting surface of the black-versus- 
reflector combination.* During the course of these 
experiments, the duct wall potential distribution was 
monitored and, as in the case of qa! there were no dis- 
cernable variations in this parameter. Accordingly, it is 
concluded that, to within the experimental uncertainty of 
gq, (46%) and EY a(t 3%) measurements, the presence of the 
radiation gage does not significantly alter the local 


thermal state of the plasma. 


* These observations are discussed further in Section A.5. 





140 


With respect to the second of the two aforementioned 
assumptions, consideration must be given to the spectral 
distribution of the incident plasma irradiation, H = H(A), 
before a suitable reflecting material, for which a, # a (A) 


Over the spectral range of H can be selected. The fol- 


i! 
lowing discussion will be restricted to the case of an 
argon plasma since the radiative wall heat flux data ob- 
Saimed in this study pertain only to this gas. 

On the basis of the theoretical work of Vakubeueee 
the emission spectrum from an atmospheric argon plasma 
is known to cover a wavelength band extending from the 
vacuum ultraviolet to the infrared. Of particular interest 
to this discussion is the vacuum ultraviolet (A < 2000 A) 
portion of the emission spectrum which results from bound- 
and free-bound transitions to the ground state (resonance 
lines and ionization continuum, respectively). As indi- 
cated by Yakubov, and on the basis of the optical depth 
and radiative recombination rate calculations performed 
by foun?’ , the radiation emitted from these transitions 
is strongly reabsorbed, and, therefore, does not contri- 
bute significantly to the balance of energy at the 
boundary of the plasma. Accordingly, it is assumed that 


there are no significant contributions to the total 


plasma radiation incident on the tube wall at wavelengths 
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oO 
less than 2000 A. This assumption is supported by the 
theoretical calculations of Maecker !1° and 1s consistent 
with the plasma radiation model employed in the non- 


L253) 


equilibrium solution by Clark as well as with the 


experimental radiation source density employed in the 


Oe te 


equilibrium solution by Bower and Incropera 
remains, however, to select a reflecting surface for which 
the monochromatic absorptance is invariant for wavelengths 


oO 
greater than 2000 A. 


A.4 Selection, Fabrication, and Optical 
Characteristics of the Reflecting Surface 


While the spectral sensitivity of a metallic sur- 
face is enhanced through polishing, the optical charac- 
teristics of polished metal surfaces are both difficult 
to reproduce and to maintain. A possible alternative 
is to adjust the optical characteristics of the surface 
by application of a thin metallic coating or thin film. 
Commonly used methods of thin film preparation include 
thermal evaporation, sputtering, and electroplating, or 
chemical deposition. The electroplating process produces 
films which are relatively porous, rough, and impure. 
This process is therefore generally ill-suited for 
applications requiring reproducible, optical quality sur- 
faces. Similarly, film deposition by sputtering is a 


relatively slow process resulting in impure and low 
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density coatings which are unsuitable for use as ultra- 
violet reflectors. Techniques employing thermal eva- 
poration in a vacuum, however, have proven to be most 
satisfactory for precise deposition of reproducible, 
optical quality, thin film coatings. The primary con- 
trol variables for this process are the rate of deposi- 
tion and the evaporator chamber pressure. 

While many metal films (e.g., aluminum, silver, 
gold, and rhodium) are highly reflecting at wavelengths 
in the visible and infrared, only aluminum films main- 
tain a high reflectance into the ultraviolet. Through 
proper preparation and aging, aluminum films will re- 
main better than 90% reflecting at wavelengths greater 
than 2000 A. For this reason, the aluminum thin fiim is 
particularly suited for the reflecting surface of the 
radiation gage. 

To achieve and to maintain the highest possible 
reflectance for the aluminum film, precise control is 
required over the following variables. 

weer time purity . 

2. Rate of deposition. 

3. Evaporator chamber pressure. 

4. Substrate temperature. 

5. Film thickness. 


6. Aging conditions. 
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The effect of each of these variables on the 
reflectance characteristics of evaporated aluminum films 


has been investigated by Madden !©8! | 


The following 
discussion provides a summary of this work. 

The purity of the aluminum has a slight effect on 
the optical characteristics of the film at ultraviolet 
wavelengths, with purer films resulting in a higher percent 
reflectance. There is little improvement, if any, in the 
percent reflectance for purities better than 99.95% and, 
for wavelengths greater than 2200 A, 99.5% pure 
aluminum films are sufficient to insure maximum reflectance. 

The condition of the oxide layer, which forms 
naturally on aluminum upon exposure to air, has a direct 
bearing on the optical properties of the aluminum coating. 
Generally, the oxide layer will grow (in air) to a thick- 
ness of from 30 to 40 A and will have little effect on 
the reflectance for wavelengths greater than 2200 A. At 
shorter wavelengths, however, the structure of the oxide 
layer is most critical, with a thinner, tighter layer 
resulting in higher reflectance characteristics for ultra- 
violet radiation. This favorable oxide layer structure 
is obtained by evaporating the aluminum film at high 
deposition rates (not less than 200 A/sec) at the lowest 


5 


possible pressure (preferably lower than 10 ~ mm Hg). 
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The substrate temperature during evaporation will 
affect the crystal size (and thus the surface roughness 
and reflectance) of the film. The larger the value of the 
substrate temperature, the lower the value of the re- 
flectance will be, particularly at the shorter wavelengths 
ee = 2200 A). However, if the substrate temperature is 
maintained below 150°C during evaporation, the reflectance 
of the aluminum film will not be affected at wavelengths 
greater than 2200 A. A similar effect results from 
increasing the film thickness beyond that required for 
maximum reflectance which, for aluminum films, is a 
thickness between 700-900 A for wavelengths greater 
than 1000 A. 

One of the more important factors influencing the 
performance of thin film aluminum reflectors is their 
age and the environmental conditions under which aging 
occurs. While there is a significant decrease in the 
percent reflectance with increasing film age at the 
shorter wavelengths, such is not the case for A > 2000 A. 
As reported by Madden, the reflectance of aluminum films 
more A > 2200 A), stored for longer than two years under 
normal laboratory conditions, remained unchanged. Regard- 


less of the age of the film, however, certain environmental 
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conditions can cause a change in its optical properties. 
imeparticular, the relative humidity and exposure of 

the film to ultraviolet radiation can produce significant 
changes in the reflectance characteristics; however, as 


reported by Hass, et a), 169] 


, these changes are minimal 
for wavelengths longer than 2000 A. 

In reviewing the previous discussion it is readily 
apparent that the reflectance characteristics of aluminum 
fms, for A > 2000 A, remain uneffected and, as shown 


by Hass, et al, 169: 


, relatively constant over a wide 
range of environmental conditions. On this basis, and 
subject to the provision that the optical characteristics 
of the film do not measureably change upon its exposure 
to the plasma, the thin film aluminum coating was 
selected for the reflecting surface of the radiation 
gage. 

Since the surface quality of a thin film is dic- 
tated by the quality of the substrate surface, special 
care was taken in the preparation of the cylindrical, 
copper surfaces prior to the deposition of the aluminum 
films. Specifically, these surfaces were ground, polished, 
and cleaned according to the following procedure. 

Grinding: 

i) 300 grit silicon carbide paper with 


water lubricant. 
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ii) 600 grit silicon carbide paper with 
water lubricant. 
Polishing: 
1) 6.0 micron diamond paste (Metadi No. 
f0-6172) with Gilelubrieant . 
1i) 0.25 micron diamond paste (Metadi No. 
10 — Gala?) (with orl lubeveant. 
1ii) 0.05 micron aluminum oxide powder with 
water lubricant. 
Cleaning: | 
1) ultrasonic cleaning with detergent agent 
(ALCONOX) . 
1i) Ultrasonic cleaning with isopropyl 


alcohol - blown dry with nitrogen gas. 


Subsequent to this procedure, the aluminum films 
were deposited with the vacuum deposition system 
(CVE-15 by Consolidated Vacuum Corporation) shown in 
Figure A.2. In an attempt to gain a uniform deposit 
over the entire semicylindrical (0.5 cm radius) surface 
of the half-segment, a special jig was designed to 
Simulate a line source of evaporation. As shown in 
Figure A.3, Enns device positions the half-segment such 
that the filament (or evaporation source) extends along 
the axis of the segment wall surface. The filament is 
placed in spring-tension so as to maintain its proper 


positon upon heating. 
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Figure A.2 Vacuum Deposition System 
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Figure A.3 
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Due to the close proximity of the substrate and 
the filament, there was some difficulty in controlling 
the thicknesses of the aluminum films. Therefore in 
order to insure that the desired 700-900 A film thick- 
nesses were achieved, a glass micro-slide, located at 
the same position from the filament as the copper sub- 
strate surface, was deposited with an aluminum film at 
the same time that each of the half-segments were film 
deposited. These glass slides were then used in a 
Twyman-Green interferometer to determine if the desired 
film thicknesses were obtained. All of the aluminum 
films were deposited at a rate of approximately 250 A/sec 
at a pressure of 2-5 x aie mn Hg. 

While the.complete spectral distribution (for 
ee 2000 A) for the monochromatic reflectance (p,) of 
the aluminum film was not measured in this study, Pp, was 
measured at several wavelengths for comparison with the 
monochromatic reflectance data reported by Drummeter and 


ae An aluminum thin film was deposited on a water- 


Hass 
cooled, copper disc (shown in Figure A.1) and this disc was 
employed for the reflectance measurements of this study. 
These measurements, which were conducted before and after 
exposure of the disc to a free argon plasma jet, were made 


: o 
over the combined wavelength range (4880 to 6328 A) of the 


argon and helium-neon laser sources employed for these 
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experiments. The results of these experiments indicate 
that there 1s no measureable change in the reflectance 
of the aluminum film after its exposure to the plasma 
environment. At each of the wavelengths investigated, 
the measured reflectance of the aluminum film was 0.98 - 
a value which is some 6-7% higher than the corresponding 
values measured by Drummeter and Hass. While the re- 
flectance measurements reported by Drummeter and Hass 
were averaged (over the wavelength interval, 2200 A < 

r < 70,000 A) and used to determine the average absorp- 
tance of the reflecting surface (01, =] - pa OO jay 
it is noted that use of the reflectance measurements of 
this study would result in values for the radiative wall 


heat flux which are within the experimental uncertain- 


ties shown in Figure IV.3. 


A.5 Presentation and Pn teisenetadicon OF 


Radiative Wall Heat Flux Measurements 

As discussed in Chapter II, it was initially felt 
that the presence of wall contamination in the arc heating 
region would preclude any meaningful radiation measurements 
in this section of the duct. For this reason, the 
radiation gage was first used in the field-free region 
(downstream from the anode), since, in this section, wall 
deposits had not been observed. With reference to 
Equations A.7 and A.9 and using values of a, = 0.98 and 


a5 = 0.07, it was determined from the calorimetric 
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measurements performed in these experiments that, of the 
total wall heat flux for the reflecting half-segment 
(qo), 2% 1s contributed by the radiative component. 

While attempts to eliminate the wall contamination 
in the arc heating region had failed, radiation experi- 
ments were nonetheless performed in this region in the 
hope that some meaningful information could be acquired. 
Upon disassembly of the duct after some four hours of high 
current (I > 100 amp) operation, it was observed that, 
while the black surface of the radiation gage and the 
copper surfaces of the adjoining segments were contamina- 
ted, the reflecting surface of the radiation gage was 
completely free of contaminants. Since the deposition 
of contaminants on the black surface precluded specifi- 
cation of its optical properties, the radiation gage 
could no longer be used in conjunction with Equation A.7 
to provide measurement of the radiative wall heat flux 
in the arc heating region. In contrast, the absence 
of contamination on the reflecting surface continues to 
permit specification of its optical properties, speci- 


Preaily oo It was therefore assumed that the radiation 


2° 
gage measurements performed in the field-free region 

could be used to delineate the radiative and convective 
wall heat flux components in the arc heating region. That 
is, of the calorimetrically measured total wall heat 


flux for the reflecting half-segment, 2% is due to 


radiation. 
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On this basis, a series of experiments were con- 
ducted with the radiation gage to determine the radiative 
and convective wall heat flux for the asymptotic region 
in argon. As indicated in the previous section, these 
experiments were further used to validate the assumption 
that the presence of the radiation gage does not signifi- 
cantly alter the local thermal state of the plasma. The 
results of this experimental investigation are presented 
in Table A.l. 

Prior to an interpretation of the data presented 
in Table A.1, it is noted that meaningful radiation 
measurements for arc currents less than 100 amp could 
not be obtained. This restriction is due to the fact 
that, at the lower arc currents in argon, the radiative 
contribution to the total wall heat flux is too small 
to permit its resolution within the estimated uncertainty 
of these measurements (+ 10%). 

The asympotic convective wall heat flux values 
shown in Table A.1 were obtained by taking 98% of the 
calorimetrically measured values for the total wall heat 
flux to the reflecting surface. The first three of the 
five separate experimental runs that were performed 
represent an attempt to determine if there exist both 
anamolies in the AT - transducer outputs and asymmetric 


conditions within the arc. While the data for the first 
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merece runs suggest that these conditions are present, the 
evidence is not conclusive since, with the exception of 
the 150 amp value for Run No. 2, all of the ar values 
agree to within the experimental uncertainty for wall 

heat transfer measurements (+ 6%). 

The remaining two experiments (Runs No. 4 and 5) 
employed a reflector-versus-reflector surface combination 
for the radiation gage, in which case the a eenewe results 
were obtained by taking 98% of the combined measured 
total wall heat flux values for both half-segments. A 
comparison between the data obtained in these two experi- 
ments indicates that the purity of the gas has an effect 
on the measured results; however, as was the case in the 
First three runs, the data agree to within the experi- 
mental uncertainties. 

In comparing the combined results for all five 
experimental runs, it is noted that, again with the 
exception of the 150 amp value for Run No. 2, the data 
agree to within + 6%. In particular, it is observed that 
such agreement exists between the data acquired with the 
black-versus-reflector surface combination and that 
acquired with the reflector-versus-reflector surface 
combination. This fact suggests that the presence of the 
radiation gage does not significantly alter the local 


thermal state of the plasma. Accordingly, the convective 
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wall heat flux data obtained with the radiation gage can 
be aeneaned with the qd measurements obtained with the 
standard copper segments in the absence of the radiation 
gage. For this purpose, the Bena values for each 

arc current condition were algebraically averaged and 

the resulting values were subtracted from the corres- 
ponding qa values (Table IV.1) to give the radiative wall 
heat flux results shown in Table A.1. The values for 

the fraction of radiative to total wall heat flux shown 


in Table A.1 were used for the comparisons in Chapters 


Gln and IV. 
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wall heat flux data obtained with the radiation gage can 
be compared with the q, measurements obtained with the 
standard copper segments in the absence of the radiation 
gage. For this purpose, the qeaaun values for each 

arc current condition were algebraically averaged and 

the resulting values were subtracted from the corres- 
ponding qs values (Table IV.1) to give the radiative wall 
heat flux results shown in Table A.1l. The values for 

the fraction of radiative to total wall heat flux shown 


in Table A.1 were used for the comparisons in Chapters 


III and IV. 
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